


CIENCIAS EXACTAS 









- YDEJA TIERRA: 


Observación, formulación y predicción 








e ANE 
TT UTOT WNVOtp cu CT 


| 


ACP TO UV 


Ino aco 


AT CL 
TIME T po LU 


FRANCISCO ODÉCIO SALES 


“ (Organizador) 





[Atena | 


Ano 2022 





CIENCIAS EXACTAS 
Y DE LA TIERRA: 


Observación, formulación y predicción 








FRANCISCO ODÉCIO SALES 
(Organizador) 


Atena 


Ano 2022 


Editora chefe 
Profº Drº Antonella Carvalho de Oliveira 
Editora executiva 
Natalia Oliveira 
Assistente editorial 
Flávia Roberta Barão 
Bibliotecária 
Janaina Ramos 
Projeto gráfico 
Camila Alves de Cremo 
Daphynny Pamplona 
Gabriel Motomu Teshima 2022 by Atena Editora 
Luiza Alves Batista Copyright O Atena Editora 
Natália Sandrini de Azevedo Copyright do texto O 2022 Os autores 
Imagens da capa Copyright da edição O 2022 Atena Editora 
iStock Direitos para esta edição cedidos à Atena 
Edição de arte Editora pelos autores. 
Luiza Alves Batista Open access publication by Atena Editora 


(DOC) Todo o conteúdo deste livro está licenciado sob uma Licença de Atribuição 
(cc) Creative Commons. | Atribuição-Não-Comercial-NãoDerivativos 4.0 
BY NC ND Internacional (CC BY-NC-ND 4.0). 


O conteúdo dos artigos e seus dados em sua forma, correção e confiabilidade são de responsabilidade 
exclusiva dos autores, inclusive não representam necessariamente a posição oficial da Atena Editora. 
Permitido o download da obra e o compartilhamento desde que sejam atribuídos créditos aos autores, 
mas sem a possibilidade de alterá-la de nenhuma forma ou utilizá-la para fins comerciais. 





Todos os manuscritos foram previamente submetidos à avaliação cega pelos pares, membros do 
Conselho Editorial desta Editora, tendo sido aprovados para a publicação com base em critérios de 
neutralidade e imparcialidade acadêmica. 


A Atena Editora é comprometida em garantir a integridade editorial em todas as etapas do processo 
de publicação, evitando plágio, dados ou resultados fraudulentos e impedindo que interesses 
financeiros comprometam os padrões éticos da publicação. Situações suspeitas de má conduta 
científica serão investigadas sob o mais alto padrão de rigor acadêmico e ético. 


Conselho Editorial 

Ciências Exatas e da Terra e Engenharias 

Prof. Dr. Adélio Alcino Sampaio Castro Machado - Universidade do Porto 

Profº Drº Alana Maria Cerqueira de Oliveira - Instituto Federal do Acre 

Profº Dr? Ana Grasielle Dionísio Corrêa - Universidade Presbiteriana Mackenzie 
Profº Dr? Ana Paula Florêncio Aires - Universidade de Trás-os-Montes e Alto Douro 
Prof. Dr. Carlos Eduardo Sanches de Andrade - Universidade Federal de Goiás 
Profº Drº Carmen Lúcia Voigt - Universidade Norte do Paraná 


[Atena 


Editora 


Ano 2022 





Prof. Dr. Cleiseano Emanuel da Silva Paniagua - Instituto Federal de Educação, Ciência e Tecnologia de Goiás 
Prof. Dr. Douglas Gonçalves da Silva - Universidade Estadual do Sudoeste da Bahia 

Prof. Dr. Eloi Rufato Junior - Universidade Tecnológica Federal do Paraná 

Prof? Dr? Érica de Melo Azevedo - Instituto Federal do Rio de Janeiro 

Prof. Dr. Fabrício Menezes Ramos - Instituto Federal do Pará 

Profº Dra. Jéssica Verger Nardeli - Universidade Estadual Paulista Júlio de Mesquita Filho 

Prof. Dr. Juliano Bitencourt Campos - Universidade do Extremo Sul Catarinense 

Prof. Dr. Juliano Carlo Rufino de Freitas - Universidade Federal de Campina Grande 

Profº Dr? Luciana do Nascimento Mendes - Instituto Federal de Educação, Ciência e Tecnologia do Rio Grande 
do Norte 

Prof. Dr. Marcelo Marques - Universidade Estadual de Maringá 

Prof. Dr. Marco Aurélio Kistemann Junior - Universidade Federal de Juiz de Fora 

Prof. Dr. Miguel Adriano Inácio - Instituto Nacional de Pesquisas Espaciais 

Profº Drº Neiva Maria de Almeida - Universidade Federal da Paraíba 

Profº Drº Natiéli Piovesan - Instituto Federal do Rio Grande do Norte 

Profº Drº? Priscila Tessmer Scaglioni - Universidade Federal de Pelotas 

Prof. Dr. Sidney Gonçalo de Lima - Universidade Federal do Piauí 

Prof. Dr. Takeshy Tachizawa - Faculdade de Campo Limpo Paulista 


[Atena 


Editora 


Ano 2022 





Ciencias exactas y de la tierra: observación, formulación y predicción 


Diagramação: 
Correção: 
Indexação: 
Revisão: 
Organizador: 


Camila Alves de Cremo 
Yaiddy Paola Martinez 
Amanda Kelly da Costa Veiga 
Os autores 

Francisco Odécio Sales 


Dados Internacionais de Catalogação na Publicação (CIP) 


C569 Ciencias exactas y de la tierra: observación, formulación y 
predicción / Organizador Francisco Odécio Sales. - 
Ponta Grossa - PR: Atena, 2022. 


Formato: PDF 


Requisitos de sistema: Adobe Acrobat Reader 
Modo de acesso: World Wide Web 


Inclui bibliografia 


ISBN 978-65-2983-908-/ 
DOI: https://doi.org/10.22533/at.ed.08/220703 


1. Ciências exatas. |. Sales, Francisco Odécio 


(Organizador). Il. Título. 





CDD 507 


Elaborado por Bibliotecária Janaina Ramos - CRB-8/9166 





Atena Editora 


Ponta Grossa - Paraná - Brasil 
Telefone: +55 (42) 3323-5493 


Wwww.atenaeditora.com.br 
contatodatenaeditora.com.br 


[Atena 


Editora 


Ano 2022 


DECLARAÇÃO DOS AUTORES 


Os autores desta obra: 1. Atestam não possuir qualquer interesse comercial que constitua um conflito 
de interesses em relação ao artigo científico publicado; 2. Declaram que participaram ativamente da 
construção dos respectivos manuscritos, preferencialmente na: a) Concepção do estudo, e/ou 
aquisição de dados, e/ou análise e interpretação de dados; b) Elaboração do artigo ou revisão com 
vistas a tornar o material intelectualmente relevante; c) Aprovação final do manuscrito para 
submissão.; 3. Certificam que os artigos científicos publicados estão completamente isentos de dados 
e/ou resultados fraudulentos; 4. Confirmam a citação e a referência correta de todos os dados e de 
interpretações de dados de outras pesquisas; 5. Reconhecem terem informado todas as fontes de 
financiamento recebidas para a consecução da pesquisa; 6. Autorizam a edição da obra, que incluem 
os registros de ficha catalográfica, ISBN, DOI e demais indexadores, projeto visual e criação de capa, 
diagramação de miolo, assim como lançamento e divulgação da mesma conforme critérios da Atena 


Editora. 


[Atena 


Editora 


Ano 2022 





DECLARAÇÃO DA EDITORA 


A Atena Editora declara, para os devidos fins de direito, que: 1. A presente publicação constitui apenas 
transferência temporária dos direitos autorais, direito sobre a publicação, inclusive não constitui 
responsabilidade solidária na criação dos manuscritos publicados, nos termos previstos na Lei sobre 
direitos autorais (Lei 9610/98), no art. 184 do Código Penal e no art. 927 do Código Civil; 2. Autoriza 
e incentiva os autores a assinarem contratos com repositórios institucionais, com fins exclusivos de 
divulgação da obra, desde que com o devido reconhecimento de autoria e edição e sem qualquer 
finalidade comercial; 3. Todos os e-book são open access, desta forma não os comercializa em seu 
site, sites parceiros, plataformas de e-commerce, ou qualquer outro meio virtual ou físico, portanto, 
está isenta de repasses de direitos autorais aos autores; 4. Todos os membros do conselho editorial 
são doutores e vinculados a instituições de ensino superior públicas, conforme recomendação da 
CAPES para obtenção do Qualis livro; 5. Não cede, comercializa ou autoriza a utilização dos nomes e 
e-mails dos autores, bem como nenhum outro dado dos mesmos, para qualquer finalidade que não o 


escopo da divulgação desta obra. 


[Atena 


Editora 


Ano 2022 





APRESENTAÇÃO 


A obra “Ciencias exactas y de la tierra: Observación, formulación y predicción” 
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Biodiversidade; Ciências Biológicas; Ciência da Computação; Engenharias; Geociências; 
Matemática/ Probabilidade e Estatística e Química. O conhecimento das mais diversas 
áreas possibilita o desenvolvimento das habilidades capazes de induzir mudanças de 
atitudes, resultando na construção de uma nova visão das relações do ser humano com 
o seu meio, e, portanto, gerando uma crescente demanda por profissionais atuantes 
nessas áreas. A ideia moderna das Ciências Exatas e da Terra refere-se a um processo de 
avanço tecnológico, formulada no sentido positivo e natural, temporalmente progressivo 
e acumulativo, segue certas regras, etapas específicas e contínuas, de suposto caráter 
universal. Como se tem visto, a ideia não é só o termo descritivo de um processo e sim um 
artefato mensurador e normalizador de pesquisas. Neste sentido, essa obra é dedicada 
aos trabalhos relacionados a pesquisa e inovação. A importância dos estudos dessa 
vertente, é notada no cerne da produção do conhecimento, tendo em vista o volume de 
artigos publicados. Nota-se também uma preocupação dos profissionais de áreas afins em 
contribuir para o desenvolvimento e disseminação do conhecimento. Os organizadores da 
Atena Editora, agradecem especialmente os autores dos diversos capítulos apresentados, 
parabenizam a dedicação e esforço de cada um, os quais viabilizaram a construção dessa 
obra no viés da temática apresentada. Por fim, desejamos que esta obra, fruto do esforço 


de muitos, seja seminal para todos que vierem a utilizá-la. 
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RESUMEN: En este trabajo, se analizó la 
capacidad de remoción de mercurio (Il) en solución 
acuosa por la biomasa del hongo Aspergillus 
niger por el método de Espectrofotometria de 
Absorción Atómica. El hongo crece hasta en 
2000 ppm del metal (20.3% de crecimiento con 
respecto al control). Se evaluó la bioadsorción 
a diferentes valores de pH (3.5, 4.5 y 5.5) y a 
diferentes tiempos. También se estudió el efecto 
de la temperatura en el intervalo de 28ºC hasta 
45ºC y la remoción a diferentes concentraciones 
iniciales de mercurio (Il) de 100 a 500 mg/L. 
La mayor bioadsorción (83.2% con 100 mg/L 
del metal y 1 g de biomasa) fue a las 24 h, a 
pH de 5.5. Con respecto a la temperatura, 
la más alta remoción fue a los 28ºC, con un 
83.2% de remoción a las 24 h, y a mayores 
concentraciones de biomasa, la remoción fue 
más eficiente (100%, 12 hy 5 g de biomasa). La 
biomasa natural mostró una excelente capacidad 
de remoción del metal in situ, 69% de remoción 
a partir de aguas contaminadas, a los 7 días de 
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incubación y 5 g de la biomasa (100 mL de agua), por lo que se puede utilizar para eliminarlo 
de aguas residuales industriales. 
PALABRAS CLAVE: Remoción, biomasa, Aspergillus niger, mercurio. 


APPLICATION OF THE BIOMASS OF THE FUNGUS Aspergillus niger FOR THE 
REMOVAL OF MERCURY IN SOLUTION 


ABSTRACT: In this work, the removal capacity of mercury (Il) in aqueous solution by 
the biomass of the Aspergillus niger fungus was analyzed by the Atomic Absorption 
Spectrophotometry method. The fungus grows up to 2000 ppm of the metal (20.3% growth 
with respect to the control). Bioadsorption was evaluated at different pH values (3.5, 4.5 and 
5.5) and at different times. The effect of temperature in the range of 28ºC to 45ºC and removal 
at different initial concentrations of mercury (Il) of 100 to 500 mg/L were also studied. The 
highest bioadsorption (83.2% with 100 mg/L of the metal and 1 g of biomass) was at 24 h, at 
a pH of 5.5. Regarding temperature, the highest removal was at 28ºC, with 83.2% removal at 
24 h, and at higher biomass concentrations, removal was more efficient (100%, 12 h and 5 
g of biomass). The natural biomass showed an excellent in situ metal removal capacity, 69% 
removal from contaminated water, at 7 days of incubation and 5 g of biomass (100 mL of 
water), so it can be used for remove it from industrial wastewater. 

KEYWORDS: Removal, biomass, Aspergillus niger, mercury. 


INTRODUCCIÓN 


El mercurio y los compuestos de mercurio están presentes en numerosos tipos 
de productos de consumo e industriales, por ejemplo: Extracción de oro y plata de las 
minas, auxiliar en la producción de químicos de cloro-álcali, en manómetros, termómetros, 
interruptores eléctricos y electrónicos, lâmparas fluorescentes, amalgamas dentales, 
aleado con otros metales, en la fabricación de pilas, como biocidas en la industria del 
papel, en pinturas o en semillas, antisépticos en productos farmacéuticos, para análisis 
químicos, y como catalizadores, para hacer más eficaz la fabricación de otras sustancias 
químicas, en pigmentos y tintes, detergentes y explosivos (Rocha-Român y cols., 2018). 
En el medio ambiente acuático, el mercurio elemental probablemente quedará unido al 
sedimento y luego será transportado por las corrientes marinas o fluviales. Una parte del 
mercurio permanece disuelta en el agua. En los sistemas acuáticos, los microorganismos 
ahí presentes de manera natural pueden transformarlo en metilmercurio, un compuesto 
organometálico que es más tóxico a dosis bajas que el mercurio en estado puro. El 
metilmercurio pasa a formar parte de la cadena alimenticia acuática, se bioacumula y se 
biomagnifica, de esta manera puede ser transportado posteriormente por las especies 
migratorias (Bell y cols., 2014). La exposición ocupacional al mercurio sigue siendo un 
problema actual para los trabajadores de muchas industrias como la mineria del mercurio, 
producción de cloro-álcali, fabricación de termómetros, lâmparas fluorescentes, baterias y 


otros productos que contienen mercurio, extracción y refinación de oro, plata, cobre y níquel 
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quienes sufren la mayor exposición son los millones de trabajadores de la extracción de oro 
artesanal y aquellos que extraen el metal en pequefa escala. Estos mineros usan mercurio 
elemental para separar el oro de los desechos de roca, generalmente en condiciones poco 
o nada controladas. Como resultado de ello, los mineros, sus familias y sus comunidades 
están muy expuestos al metal, con sus muy graves consecuencias para su salud (Bell y 
cols., 2014). 


LA PRODUCCIÓN DE MERCURIO EN MÉXICO 


Los minerales de mercurio se encuentran localizados principalmente en las 
siguientes entidades federativas (en particular en el noroeste y centroeste del territorio): 
Aguascalientes, Chihuahua, Coahuila, Durango, Guanajuato, Hidalgo, México, Morelos, 
Nuevo León, Querétaro, San Luis Potosí, Sinaloa, Sonora, Tamaulipas, Tlaxcala y 
Zacatecas, ubicândose los principales yacimientos en los estados de: San Luis Potosí, 
Zacatecas, Querétaro, Guanajuato y Guerrero (Figura 1, cuadros negros) (Consejo de 
Recursos Minerales 2002). 





Eu | Mercurio | 


Figura 1.- Modificada de: Principales yacimientos de mercurio en la República Mexicana. (Consejo de 
Recursos Minerales 2002). 





Se tienen registros de la producción de mercurio en México desde 1891, y en los 
Últimos ahos se ha mostrado la misma tendencia descendente que la producción mundial. 
Entre 1920-1929 su producción fue mínima y alcanzó su máximo de 1,118 toneladas en 
1942; en 1991 se produjeron 340 toneladas de mercurio, mientras que en 1994 sólo 11 
toneladas. A partir de 1995 no se tienen registros de su extracción minera (Consejo de 
Recursos Minerales 2002). Sin embargo, la producción secundaria de mercurio continua, 
en particular, al recuperarlo como subproducto de la extracción de plata y oro a partir de 
jales (deshechos) mineros en el estado de Zacatecas, donde se producen alrededor de 20 


toneladas anuales. Entre las plantas productoras de minerales metálicos registradas en 
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la Câmara Minera de México, se han descrito cuatro, que se dedican a la producción de 
mercurio a partir del beneficio de minerales y de jales de mercurio, las cuales se encuentran 
ubicadas en los estados de San Luis Potosi, Durango y Zacatecas. Asimismo, se tiene 
conocimiento que hasta hace algunos afios operaban dos plantas más en el estado de 
Zacatecas, que benefician jales por el proceso de lixiviación, utilizando hiposulfito de 
sodio, para obtener oro y plata como producto y mercurio como subproducto, y debido a la 
reducción de los usos comerciales del metal tanto a nivel internacional como en México, y 
por la saturación del mercado nacional, se ha reducido su producción hasta haber dejado 
de producirlo desde 1995 a la fecha (INE 2001). 

Por lo anterior, la situación que se vive en México con respecto a la contaminación 
y/o producción de mercurio, es contradictoria, pues a pesar de las disposiciones normativas 
formuladas, que limitan las emisiones de mercurio a la atmósfera y al agua y controlan la 
eliminación de residuos que contienen mercurio, el elemento no ha sido regulado como 
producto comerciable y se ha hecho poco por fomentar entre la población una adecuada 
conciencia en cuanto a la exposición al mercurio y la reducción de sus riesgos, y de acuerdo 
con el Artículo 7º de la Ley General para la Prevención y Gestión Integral de los Residuos, 
los residuos peligrosos son de competencia de la federación, mientras que la NOM- 
052-SEMARNAT-2005 establece las características, el procedimiento de identificación, 
clasificación y los listados de los residuos peligrosos, incluyendo aquellos que contienen 
mercurio, incluyendo diversos residuos con mercurio, clasificados como residuos peligrosos, 
y se establece como criterio de toxicidad la concentración de diversas sustancias en el 
extracto PECT (cuyo proceso se establece en la NOM-053-SEMARNAT-1993), siendo de 
0.2 mg/L para el caso de mercurio. 


ANTECEDENTES 


La bioadsorción es un proceso fisicoquímico que incluye los fenómenos de adsorción 
y absorción de moléculas e iones. Este método tiene como objetivo principal la remoción 
de metales pesados en aguas residuales prevenientes del sector industrial, usando como 
sorbentes diferentes materiales de origen biológico (vivos o muertos), tales como: algas, 
hongos, bacterias, cáscaras de frutas, productos agrícolas y algunos tipos de biopolímeros. 
Estos materiales son de bajo costo y se encuentran en gran abundancia en la naturaleza, 
además, su transformación a biosorbente no es un proceso costoso (Acosta Rodríguez 
y cols., 2012; Tejada Tovar y cols., 2015). El proceso de bioadsorción involucra una fase 
sólida (biomasa) y una fase líquida (agua) que contiene disueltos la sustancia de interés 
que será adsorbida (en este caso, los iones de los metales pesados). 

Investigaciones recientes sobre la adsorción de iones de metales pesados en 
aguas residuales, realizadas a nivel mundial, revelan la capacidad de adsorción de 
diferentes biomasas naturales como Aspergillus versicolor (Das et al., 2007); carbones 
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activados (Colpas et al., 2016); Omphalina sp. (Purwaningsih et al., 2016); Escherichia 
coli, Pseudomonas putida, Serratia marcescens y Xanthomonas sp. (Ashraf y Essa, 2012); 
Phanerochaete chrysosporium (Bashardoost et al., 2010); Rhizopus oligosporus (Ozsoy 
2010); silica biogénica modificada con L- cisteíina (Chaves et al., 2011); diferentes hongos 
contaminantes ambientales (Martinez Juárez et al., 2012); y el hongo Mucor hiemalis 
(Hoque y Fritscher, 2016). Por lo anterior, es de gran interés analizar otros microorganismos 
capaces de eliminar el mercurio (Il) de aguas y suelos contaminados con dicho metal, 
como la biomasa del hongo contaminante ambiental Aspergillus niger, microorganismo que 
podria remover el mercurio por interacción con los componentes estructurales de la pared 
celular, como los grupos carboxilos y amino. Además, una ventaja del uso de este hongo 


como biomasa, es que se puede propagar de manera muy fácil y económica. 


MATERIAL Y MÉTODOS 


Bioadsorbente utilizado y obtención de la biomasa fúngica 


Se trabajó con una cepa de 4. niger que crece en 2000 mg/L de mercurio (Il), aislada 
del aire contaminado de una zona adyacente a la Facultad de Ciencias Químicas de la 
UASLP, San Luis Potosí, S.L.P., México, en el mes de junio de 2015; y se identificó en 
función de sus características morfológicas, tanto macroscópicas como microscópicas 
(López Martínez et al., 2004). La cepa fúngica se mantuvo de manera rutinaria en agar 
de dextrosa y papa. Para la obtención de la biomasa, se sembraron 1 x 10º esporas/mL 
del hongo en matraces Erlenmeyer de 1 L, que contenían 600 mL de caldo tioglicolato, 
incubândose durante 5 días a 28 “Cy 100 rpm. Después, la biomasa se obtuvo por filtración 
y se lavó con 200 mL de agua tridesionizada. El producto se molió en una licuadora y se 
almacenó en frascos âmbar hasta su uso.Soluciones de Mercurio (Il). 

Se trabajó con 100 mL de una solución de 100 mg/L de mercurio (Il) obtenida por 
dilución de una solución patrón de 1 g/L preparada en agua tridesionizada a partir de cloruro 
de mercurio. Se ajustó el pH de la dilución a analizar con H,SO,1 Mo NaOH 1 M, antes de 
adicionarlo a la biomasa. 

Prueba de Resistencia al Mercurio 


Se inocularon 1 x 10º esporas/mL en 12 matraces Erlenmeyer de 250 mL, conteniendo 
100 mL de medio mínimo de Lee modificado (LMM), conteniendo 1.0 % de agar, 0.25 % 
de KH,PO,, 0.20 % de MgSO,, 0.50 % de (NH ),SO,, 0.50 % de NaCl y 0.25 % de glucosa 
como fuente de carbono, y complementado con concentraciones crecientes de mercurio 
(II), como cloruro de mercurio (0-2000 mg/L), incubando a 28ºC a 100 rpm durante 7 días. 
Después, se cosechó el sobrenadante en un tubo graduado, previamente pesado, y se 
centrífugo a 3000 rpm/O min, desechando el sobrenadante. El paquete celular se secó 
a 80ºC, durante 12 h, y se pesó el tubo, determinando por diferencia el peso seco de la 
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muestra, comparando el crecimiento con un control crecido en las mismas condiciones 
sin la adición de mercurio (Il). Todos los experimentos se realizaron mínimo 3 veces por 


duplicado. 


Estudios de remoción de mercurio (II) 


La concentración de mercurio (Il) en solución, se determinó por Espectrofotometría 


de Absorción Atómica por generación de hidruros (SSA, 1994). 


RESULTADOS 


En la Figura 2 se muestra que, a medida que se aumenta la concentración del 
metal (200-2000 mg/L) en la solución, el crecimiento de A. niger en el medio de Lee, 
disminuye en relación con el control, pues con 2 g/L de mercurio (Il), el crecimiento fue de 
aproximadamente un 20.3 % de crecimiento con respecto al control (85,4 mg de peso en 
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Figura No. 2. Crecimiento en peso seco de A. niger con diferentes concentraciones de mercurio (Il), 
1X10º esporas/mL, 28ºC, 7 días de incubación y 100 rpm. 


Con respecto al efecto del pH y tiempo de incubación sobre la remoción de 100 mg/L 


de mercurio (Il) por la biomasa del hongo A. niger, se observó que la mayor bioadsorción 
ocurre a las 24 horas y a un pH de 5.5, con una remoción del 83.2 % (Figura 3). 
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Figura No. 3. Efecto del tiempo de incubación y pH sobre la eliminación de mercurio (Il) por la biomasa 
de A. niger 100 mg / L Hg (Il), 100 rpm, 28ºC, 1 g de biomasa. 


Con relación al efecto de la temperatura de incubación, la temperatura óptima para 


la remoción de mercurio (Il) fue de 28 “C, con un porcentaje de remoción de 85 %. Dichos 
resultados se ilustran en la Figura 4. 
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Figura No. 4. Efecto de la temperatura de incubación sobre la remoción de mercurio (Il) por la biomasa 
de A. niger 100 mg / L de Hg (Il), 100 rpm, 28ºC. 1 g de biomasa. 


Además, la concentración inicial del metal influye en su remoción, pues el porcentaje 
de remoción es menor a medida que se aumenta la concentración del contaminante, debido 
a que, al haber mayor cantidad del metal disponible, hay menor interacción de las moléculas 
del contaminante con los sitios de adsorción de la biomasa (Figura 5). 
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Figura No. 5. Efecto de la concentración inicial de metal sobre la remoción de mercurio (Il) por la 
biomasa de A. niger. 100 rpm, 28º C. 1 g de biomasa. 


Por otro lado, al aumentar la concentración de la biomasa, aumenta la eficiencia 
de remoción, pues con 5 g de biomasa, el tiempo de remoción es menor para el metal 
analizado (Figura 6). 
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Figura No. 6. Efecto de la concentración de biomasa sobre la remoción de 100 mg/L de mercurio (Il) 
por la biomasa de 4. niger. 100 rpm, 28ºC. 


Finalmente, se evaluó la remoción del metal durante un periodo de incubación de 7 
dias en aguas contaminadas naturalmente (concentración inicial 183 mg/L), obteniendo una 
remoción del 69 % al final del experimento (Figura 7). Las muestras fueron obtenidas del 
Tanque Tenorio, localizado al este de la ciudad de San Luis Potosí, México. Originalmente, 
era una depresión natural y desde finales de los aos setenta del siglo XX, se empezó a 
utilizar de manera continua para descargas industriales y en menor grado, para aguas 
residuales domésticas (López Álvarez, 2012). Esas aguas se han utilizado para irrigación 
de cultivos en zonas aledafias, impactando en gran medida la salud y los acuíferos 
subterrâneos de la zona (COTAS, 2005). Se identificó que esas fuentes de contaminación 
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contienen metales pesados como: Hg, Ba, Sr, Cd, Pb, Ag, Rb, Co, Cu, Fe and As, y algunos 
de ellos exceden los límites permisibles para el consumo humano (NOM-127- SSA1-1994). 
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Figura No. 7. Eliminación de mercurio (Il) en desechos industriales incubados con 5 g de biomasa. 
28ºC, 100 rpm, pH 8.2, 100 mL de agua contaminada, (183 mg/L). 


DISCUSIÓN 


La cepa fúngica creció un 20.3% con 2 g /L de mercurio (Il), con respecto al control 
sin metal. Diferentes microorganismos que son resistentes a este metal han sido aislados 
de diferentes sitios contaminados, Giovanella y cols., (2017), reportaron que las bacterias: 
Pseudomonas putida, Pseudomonas sp., Alcaligenes faecalis y Brevundimonas sp., que 
se aislaron de residuos y efluentes contaminados con diferentes metales, crecieron con 
10 uM de mercurio. Mahbub y cols., (2016), aislaron Sphingobium, una cepa bacteriana 
altamente resistente al mercurio de suelo contaminado, con valores estimados de CE5O 
de 4.5 mg/L y 44.15 mg/L y valores de MIC de 5.1 mg/L y 48.48 mg/L; Koushalshahi y 
cols., (2012), aislaron la bacteria Streptomyces sp., a partir de sedimentos marinos en 
diferentes regiones del Mar Caspio, que crece a 20 y 40 mg/L de mercurio. Sin embargo, el 
mecanismo de tolerancia en esta cepa de hongo de A. niger no se conoce. 

Uno de los problemas ambientales más importantes que se relacionan con 
la contaminación de agua y suelo en el país, es el vertido directo de metales pesados 
altamente contaminantes, lo cual constituye un peligro para el equilibrio y el funcionamiento 
natural de los ecosistemas acuáticos. Uno de estos metales es el mercurio, utilizado en 
procesos de extracción y purificación de oro artesanal. La bioadsorción, a través del uso 
de microorganismos y residuos agroindustriales, ha demostrado ser una alternativa para 
el tratamiento de matrices acuáticas contaminadas por metales, es una tecnologia limpia, 
eficiente y económica para los procesos de remediación ambiental. Los principales factores 
asociados al proceso de adsorción de mercurio (Il) por la biomasa del hongo 4, niger, 
establecieron que este biomaterial tiene un alto potencial en la remoción de este metal, 
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logrando una remoción máxima del 83.2% (100 mg/L y 1 g de biomasa) a las 24 h, pH 
de 5.5 y 28ºC. La evidencia experimental muestra un fuerte efecto de las condiciones 
experimentales. Los valores máximos de la capacidad de bioadsorción mostraron que este 
biosorbente, es muy efectivo en la recuperación y/o eliminación de iones de mercurio de 
sistemas acuáticos. 

Finalmente, se observó que, a los 7 dias de incubación con la biomasa fúngica, 
la concentración del metal de las muestras de aguas residuales disminuye un 69%. En 
la literatura se ha reportado la capacidad del quitosano impregnado con bromuro de 
etilhexadecildimetilamonio para la eliminación de cadmio y mercurio, preparando dos 
efluentes sintéticos diferentes que conteniían Cd (Il) en la primera, y Hg (Il) en la segunda 
(composición similar a la de los efluentes reales) y se determinó la eficiencia de adsorción 
para estos iones metálicos. Los resultados muestran una eliminación de 86.4% y 73.8% 
para Cd (ll) y Hg (ll), respectivamente (Shekhawat y cols., 2017). La utilización de 
nanopartículas de selenio biogénicas y goetita elimina un 81.2% del mercurio, a partir de 
agua subterránea contaminada (Wang y cols., 2018). La remoción del 99% del metal de 
soluciones acuosas, después de 10-48 h posteriores a la exposición inicial a mercurio (Il), 
del agua subterrânea de Munich, con una nueva cepa Mucor hiemalis EH8 que acumula 
mercurio (Hoque y Fritscher, 2016). Finalmente, en la Tabla No.1, se muestra la capacidad 


de remoción de mercurio (Il) por diferentes biomasas. 


Adsorbente pH Capacidad de adsorción Referencia 
mg/g 

Aspergillus niger ha 3.2 Karunasagar, y cols., 2003 

Aspergillus versicolor 5.0-6.0 75.6 Das y cols., 2007 

Different fungal 5. 78.3-95.3 Martinez-Juárez y cols., 2012 

biomasses 

Eucaliptus leaves 5.0-6.0 129.87 Alimahammadi y cols., 2017 

Impregnated chitosan 3.0 43.3 Shekawat y cols., 2017 

Fhanerochaete 6.0 72.46 Bashardoost y cols., 2010 

chrysosporium 

Pleurotus sapidus 6.0 207.89 Yalcinkaya y cols., 2002 

Rhizopus oligosporus 60 Sds Ozsoy, 2010 

Silica gel modified 3.0-7.0 0.20 Chaves y cols., 2011 

Aspergillus niger 2.5 83.2 Este trabajo 


Tabla No 1. Comparación de la capacidad de adsorción de Hg (Il) con diferentes bioadsorbentes. 


Como se observa en la tabla anterior, la cepa utilizada en este trabajo tiene buena 
eficiencia de remoción, y es de la misma especie utilizada por Karunasagar, y cols., 
(2003), y es mucho más eficiente en la remoción del metal, lo cual puede deberse a que 
ellos utilizaron un bioadsorbente preparado a partir de la biomasa del hongo tratada con 
hidróxido de potasio al 5% y caliente, y el metal utilizado fue inorgânico (HG?*) y metil 


mercurio (HgC!º). 
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Además, algunos hongos que tienen la capacidad de remover mercurio en solución, 
también pueden eliminar otros metales como el cromo (VI) y el arsénico (Ill y V) (Acosta- 
Rodríguez y cols., 2017; Santos y cols., 2017), y estos resultados confirman la capacidad 
de esta biomasa para ser utilizada en la biorremediación de metales pesados de nichos 
contaminados, como la laguna del Tanque Tenorio descrita anteriormente. 

Por otro lado, el pH y la temperatura son parâmetros críticos en la remoción de 
metales pesados, en este trabajo se encontró que el pH y la temperatura óptimos para la 
remoción de mercurio (Il) fuero 5.5 y 28ºC, con un porcentaje de remoción de 85%, y con 
respecto a su remoción en presencia de otros metales pesados, éstos si pueden influir en la 
remoción del mercurio, aunque no de manera significativa, sobre todo el del pH, pues otros 
metales pesados se remueven a otros pH's, como la temperatura es parecida a la media 
ambiental, ésta no influiria de manera significativa en la remoción. 

Finalmente, una limitante para su aplicación en lugares contaminados, seria el 
tiempo de contacto de la biomasa con el metal, porqué se requieren 24 horas de contacto 
para tener una buena eliminación del mismo, y en suelos puede influir el viento, la Iluvia y la 
presencia de basura y animales, y en aguas, la mayoría de éstas son aguas corrientes los 
cual impediria un tiempo de contacto prolongado con la biomasa, y por lo tanto la eficiencia 
de remoción puede disminuir, aunque nosotros pensamos que se puede aumentar la 
cantidad de biomasa para tratar de optimizar el proceso, o inicialmente trabajar en aguas 


estancadas. 


CONCLUSIONES 


1. La remoción de Hg (Il) por la biomasa de A. niger fue de 83.2%, a pH 5.5, 28ºC, 
24 h, 100 rpm y 100 mg/L. 


2. La biomasa del hongo 4. niger demostró una excelente capacidad de remoción 
del metal in situ, con un 68% de remoción en aguas industriales contaminadas, a los 
7 días de incubación, 5 g de biomasa y con 100 mL de agua, a un pH de 8.2 


3. Las condiciones óptimas para la bioadsorción del metal con la biomasa de 4. 
niger fueron: T = 28ºC y pH de 5.5, 100 rpm, 24 h de incubación, 100 mg/L del metal 
y 5 g de biomasa, mientras que en el agua contaminada fueron T= 280C, pH 8.2, 100 
rpm, 7 días de incubación, 183 mg/L del metal y 5 g de biomasa. 


4. El proyecto es viable para aplicarse en sitios contaminados, ya que se comprobó 
con agua contaminada naturalmente, obteniendo una buena adsorción del metal, 
utilizando las condiciones descritas en este trabajo, y de preferencia en aguas 
estancadas. 
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ABSTRACT: In the present research work, a 
study was carried out on the quality of water for 
human consumption, in the central region of the 
P'urhépecha plateau of the state of Michoacán, 
Mexico, monitoring the variables of hardness, 
TDS, temperature, pH and concentration. of As 
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MESETA P'URHEPECHA 


and Pb in different sampling sites and in different 
seasons of the year, which were compared with 
parameters established by the Official Mexican 
Standard  NOM-127-SSA1-1994,  NOM-250- 
SSA1-2014, Environmental Protection Agency of 
the United States of America (USEPA) and those 
established by the World Health Organization 
(WHO). The results obtained show a pH that is 7 
+ 0.5, the temperature oscillates between 18 and 
22º €, the hardness between 60 and 200 ppm, 
unlike a site that shows 600 ppm. SDT between 
40 and 220 ppm and a sample site of 400 ppm. 
On the other hand, the Pb concentration range 
was found between 0.15 and 0.5 ppm, while As 
concentrations ranged between 0.01 and 0.1 
ppm, depending on the study site. 

KEYWORDS: Water quality, hardness, TDS, 
lead, arsenic. 


WATER QUALITY IN THE CENTRAL 
REGION OF THE MESETA 
P'URHÉPECHA 


RESUMEN: En el presente trabajo de 
investigación se realizó un estudio sobre la 
calidad del agua para consumo humano, en la 
región centro de la meseta P'urhépecha del 
estado de Michoacán, México, monitoreando 
las variables de dureza, SDT, temperatura, pH y 
la concentración de As y Pb en diferentes sitios 
de muestreo y en diferentes estaciones del afo, 
los cuales fueron comparados con parâmetros 
establecidos por la Norma Oficial Mexicana 
NOM-127-SSA1-1994,  NOM-250-SSA1-2014, 
Agencia de Protección Ambiental de los Estados 
Unidos de América (USEPA) y los establecidos 


Capítulo 2 


por la Organización Mundial de la Salud (OMS). Los resultados obtenidos muestran un pH 
que se encuentra en 7 + 0.5, la temperatura oscila entre los 18 y 22 ºC, la dureza entre los 
60 y 200 ppm a diferencia de un sítio que muestra 600 ppm. SDT entre 40 y 220 ppm y un 
sitio de muestreo de 400 ppm. Por otra parte, el rango de concentración de Pb se encontró 
entre 0.15 y 0.5 ppm, mientras que las concentraciones de As oscilan entre 0.01 y 0.1 ppm, 
dependiendo del sitio de estudio. 

PALABRAS CLAVE: Calidad del agua, dureza, SDT, plomo, arsénico. 


INTRODUCCIÓN 


Hablar de la calidad del agua para el consumo humano hoy en día es un tema 
de crucial importancia, debido que el vital líquido es un medio donde se Ilevan a acabo 
múltiples reacciones bioquímicas. Para el caso de las reacciones metabólicas en el ser 
humano, ésta permite el transporte de sustancias en forma iónica dentro de las células y 
entre ellas, lo que permite ser un medio indispensable para la vida humana (Voet, Voet & 
Pratt, 2008). Para el caso de la agricultura su utilidad y aplicación resulta ser básica y en 
esta misma vertiente para la industria sus aplicaciones se diversifican en múltiples usos, 
dependiendo del giro de las mismas. Sin embargo, dada su utilidad de este compuesto vital 
se ha observado que cada vez se ve una mayor contaminación, debido a la combinación de 
diferentes especies químicas, por ejemplo: nitratos, cloruros, oloruros, elementos pesados 
entre varios otros, por ello es importante monitorear aspectos como la potabilidad, la 
presencia de materia orgânica daíina, la alcalinidad, cantidad excesiva de solidos disueltos, 
dureza, conductividad entre otras variables, (Peréz-López, E., 2016). 

La calidad del agua por otro lado, implica parâmetros permisibles que no dafen la 
salud en el cuerpo humano. Particularmente el agua para beber, es mayormente consumida 
si ésta pasa por un proceso de purificación. Una clasificación del agua para consumo 
humano es la siguiente: agua potable cruda y agua de mesa, donde la primera es agua 
clorada y con excedentes de sales, por otro lado, el agua de mesa, ha recibido un proceso 
de purificación a través de filtros de sedimentos, ósmosis inversa o cualquier otro proceso 
(Vidales-Olivos, 2000), el principal objetivos de dichos procesos es disminuir la cantidad de 
solidos disueltos, dureza y algunos otros parâmetros establecidos por normas nacionales 
NOM-250-SSA1-2014 e internacionales. Sin embargo, existen regiones de la república 
mexicana donde la población hace uso del vital líquido para beber directamente de 
manantiales y pozos profundos sin someterla a ningún proceso de purificación ni cloración. 

Por otro lado, el agua que proviene de manantiales y pozos profundos puede 
contener algunos elementos pesados, como Arsénico, Plomo, Mercurio, Cadmio y varios 
otros elementos en forma de iones y sus diferentes compuestos, esto depende del subsuelo, 
minerales de contacto e incluso por contaminación de actividades por el ser humano, lo que 
implica un riesgo para la salud humana (efsa, 2015; Huang et. al, 2014). 

En la región centro de la meseta P'urhépecha, ubicado en el estado de Michoacân, 
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México, se encuentran localizados los municipios de Nahuatzen, Cherán y Paracho, donde 
la población local extrae el agua de manera cotidiana de pozos profundos y manantiales 
consumiendo sin previa potabilización y/o purificación, por ello se vuelve importante la 
evaluación de la calidad del agua en fuentes de abastecimiento de los municipios de interés 
y de estudio de la región centro de la meseta P'urhépecha. 

En esta investigación el objetivo principal es evaluar y determinar la calidad del agua 
en manantiales y pozos profundos ubicados en la región centro de la meseta P'urhépecha a 
través del análisis fisicoquímico y la determinación de metales pesados. Las variables para 
este estudio fueron: dureza, cantidad de solidos disueltos, pH, temperatura, concentración 
de As (arsénico) y la concentración de Pb (plomo). 


MÉTODOS 


Para llevar a cabo esta investigación y el cumplimiento de los objetivos se 
establecieron los siguientes pasos: 


a) Toma de muestras de agua de manantiales y pozo profundo. 


Se identificaron doce sitios de muestreo; tres en el municipio de Nahuatzen, seis en 
el municipio de Cherán y tres en el municipio de Paracho, Se tomaron dos muestras por 
sitio en viales de 50 ml en las diferentes estaciones del aho. 


b) Análisis fisicoquímico; dureza, cantidad de solidos disueltos, pH y temperatura. 


Por cada muestra tomada en in situ se midió la dureza, solidos disueltos, pH y 
temperatura con equipos portátiles (electrodo de dureza, de TSD, pH-chimetro y termómetro 
laser). 

c) Análisis fisicoquímico por emisión atómica por plasma para determinar metales 
pesados. 

Las muestras tomadas fueron analizadas en un espectrómetro de emisión por 
plasmas modelo 4100 MP-AES marca Agillen. 

d) Análisis estadístico de las zonas muestreadas. 
Una vez analizadas las muestras en sus diferentes variables por cada estación del 


ano, se promedió el resultado y posteriormente se obtuvieron graficas en donde se aprecia 
de forma comparativa cada resultado de los manantiales y pozos profundos estudiados. 


RESULTADOS 


A continuación, se muestra una serie de gráficos en la que se puede apreciar el 
monitoreo y evolución de la cantidad del agua medida en ppm de las siguientes variables: 
dureza, cantidad de solidos disueltos, pH, temperatura, concentración de As y la 
concentración de Pb. Cabe resaltar que estos resultados dan información relevante de la 
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calidad del agua de los diferentes sitios que fueron muestreados. 
pH 

En la figura 1, se puede apreciar el pH promedio del agua de cada uno de los sitios 
donde fue muestreada en in situ. La mayoría de los pozos y manantiales se encuentra a 
un pH de 7 más, menos 0.5 unidades de diferencia entre algunos sitios. Este resultado se 
encuentra acorde con la norma permisible mexicana NOM-127-SSA1-1994 y la NOM-250- 
SSA1-2014, el cual establecen que el agua potable debe estar en un rango de pH entre 
6.5-8.5. Cabe mencionar que la toma de muestras se realizó en diferentes temporadas del 


ano y la variación que se encontró fue mínima (entre 0.4 y 0.5). 


O—O 





Figura 1. pH promedio del agua de zonas muestreadas. 


Temperatura 


Se encontró que la temperatura promedio varia según el sitio de muestreo. 
De acuerdo a la gráfica de la figura 2, los sitios de San Juan (P), Pozo 4 2 y Pozo & 4 
pertenecientes a la localidad de Paracho diferenciado por la letra (P), el cual, mostraron 
un rango de temperatura entre los 19 y 22 ºC, por otra parte, para los sitios de muestreo 
de la localidad de Cherán (diferenciado con la letra C): Manantial Cofradia (C), Manantial 
Tendeparacua (C), Pozo Plan (C), Pozo del Rebombeo (C), y Manantial Cumitzaro (C), la 
temperatura promedio es muy similar a la localidad de Paracho, es decir; que esta, entre 
los 18y 21 ºC, a diferencia del Pozo del Barrio Tercero (C) donde el promedio se encontró 
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entre los 25 ºC, este mismo resultado fue similar al que se obtuvo con los sitios de la 
comunidad de Nahuatzen (diferenciados con la letra N), de los tres puntos de muestreo; 
San Francisco Cerro Capen (N), Pozo Profundo (N) y Ojo de agua del Plan. 


Temperatura (C) 





Figura 2. Temperatura promedio del agua de zonas muestreadas. 


Dureza 


Cuando se habla de dureza se hace mención a la cantidad de carbonato de calcio 
(CaCO,) disuelto en el agua en mayor concentración y en menor el carbonato de magnesio 
(MgCO.). En este sentido el agua puede tener diferentes connotaciones, agua dura, agua 
parcialmente dura o blanda. 

La Organización Mundial de la Salud (OMS), define la clasificación como agua blanda 
la que presenta concentraciones inferiores a 60 mg/L de Carbono de Calcio (CalCO.), entre 
61 y 120 mg/L, moderadamente dura, entre 121 y 180 mg/L dura y aquella con valores 
superiores a 180 mg/L muy duras (Instituto Nacional de Epidemiológica y Microbiológica, 
1992). 

En este estudio en particular se tomó, en cuenta esta variable, puesto que es un 
parâmetro fisicoquímico de importancia para la calidad del agua de consumo humano 
(potable). Para lo cual se encontraron los siguientes resultados. En los pozos de la localidad 
de Cheran la dureza, figura 3, muestra que se encuentran entre 80 y 100 ppm (o mg/L), 
para algunos sitios, mientras que para otro oscila entre los 100 y 150 ppm (manantial 
Tendeparacua), a diferencia del sitio conocido como Manantial Cumitzaro donde la 
concentración promedio de los carbonatos esta entre los 600 ppm. 
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Para las zonas de muestreo de Nahuatzen muestran un rango promedio entre los 60 
y 100 ppm, La mayoría de los sitios de muestreo muestran que los resultados se encuentran 
dentro de la norma permisible establecidos por la NOM-250-SSA1-2014, el cual establece 
que la dureza debe estar a 500 ppm como límite máximo permisible. En contraste al sitio 
de muestreo del Manantial Cumitzaro del municipio de Cherán, donde arrojo valores por 
arriba de la norma, 600 ppm aproximadamente, esto indica un foco de atención para los 
habitantes que hacen uso de tal recurso hídrico. Por otro lado, en el municipio de Paracho 
en rango de los sitios de muestreados escila entre los150 y 200 ppm. 


Ojo de agua del pilon (NJ) 
Pozo Profundo (NJ/77777) 


San Francisco cerro Capen 7 


Manantial Cumitzaro (7/1/1111 
Pozo del barrio tercero (C77777 
Pozo de rebombeo (077777 
Pozo del plan (C)777777 
Manantial Tendeparacua (C)7777777 
Manantial Cofradia (C)7777777) 
Pozo 4 (777/77 
Pozo É 2. (P777777) 
San Juan (P MU 


DU DU JU 4UU DU 6OU 
Dureza (ppm) 


Figura 3. Dureza promedio del agua de zonas muestreadas. 


Por otro lado, la aceptación de la dureza varía dependiendo de la aceptación del 
público, incluso cuando se habla de agua embotellada para beber, esta depende de las 
especificaciones de calidad al cual están sometidas. Normalmente el rango de sabor del ión 
calcio esta entre los 100 a 300 ppm y el de magnesio es probablemente menor al de calcio. 
Segun GEMS (1992) si la cantidad de magnesio es superior a 125 ppm en el agua de beber 
esta puede tener efectos laxantes. En este contexto si el exceso de carbonato de calcio en 
agua de beber está presente, esta puede generar cálculos renales en el organismo humano 
(Mora y Chamizo, 2007). Así mismo Melitón y Rosado (2008) concluyen que el consumo de 
agua con niveles por arriba de 120 ppm de carbonato de calcio son factores importantes 
en la producción de cálculos renales en las vias urinarias. Por otro lado, investigaciones 
realizadas mencionan que las aguas blandas y duras han sido relacionadas con problemas 
cardiacos (OMS, 1983) y cálculos renales (Evans y Taal, 2011; Mora y Alfaro, 1999). En 
acuerdo con los autores antes mencionados se muestran evidencia de la importancia que 


la dureza en agua potable es una variable que debe ser monitoreada de manera constante 
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y más aún si el agua es de meza. 
Sólidos Disueltos Totales (SDT) 


Los sólidos son partículas de compuestos que se encuentran diluídas en agua natural, 
se caracterizan por permanecer en estado sólido después de la evaporación del líquido. 
Básicamente existen de dos tipos sólidos suspendidos totales (SST) y los sólidos disueltos 
totales (SDT), correspondientes a residuos no filtrables y filtrables, respectivamente. 

En la figura 4, se puede apreciar los resultados obtenidos de los sitios de muestreo, 
donde se muestra la cantidad promedio de SDT, para los sitios de la comunidad de 
Nahuatzen se observa que están en un rango de entre 60 y 90 ppm, mientras que para 
los de la comunidad de Cheran el rango se encuentra entre los 50 y 70 ppm a diferencia 
del sitio Manantial de Cumitzaro donde se muestra una cantidad promedio de 400 ppm. 
Por otro lado, en la localidad de Paracho, el rango de SDT está entre los 140 y 200 ppm. 
Cabe mencionar que estos resultados muestran evidencia que se encuentran dentro de la 


cantidad máxima permisible. 


Ojo de agua del pilon (N)JKWN 
Pozo Profiundo (NJ 
San Francisco cerro Capen (N)JAANN 
Marat Câm (SSIS 
Pozo del barrio tercero (CNN 
Pozo de rebombeo (ON 
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Manantial Tendeparacua (ON 
Manantial Cofradia (CAN 
Pozo & 4 (P AIN 
Pago ff 2 (PANA 
San Juan (P ANN NES 
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Figura 4. Promedio de sólidos disueltos en agua de las zonas muestreadas. 


En contraste, la NOM-127-SSA1-1994 establece una máximo permisible de SDT de 
1000 ppm (mg/L), mientras que la agencia de protección ambiental de los Estados Unidos 
de América (USEPA), sugiere un máximo de 500 ppm en agua potable, esto debido a que 
si la cantidad de SDT está arriba de lo sugerido, el agua puede presentar apariencia turbia 
y con sabor distinto, incluso las personas no acostumbradas al agua con alto contenido de 
SDT pueden experimentar irritación gastrointestinal al beber agua con esta característica 
(Sigler y Bauder, 2017). 
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Pb (Plomo) 


En la siguiente figura 5, se puede apreciar los resultados obtenidos de la 
concentración de Pb encontradas en los diferentes manantiales y pozos profundos de las 
localidades antes mencionadas. 
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Figura 5. Concentración promedio en ppm (mg/L) de Pb en sitios de muestreo. 


Los pozos profundos y manantiales de lo localidad de Paracho muestran 
concentraciones promedio de 0.46 y 0.5 ppm (mg/L), en los de Chéran se encuentran en 
un rango de 0.25 y 0.48 ppm, para los de la localidad de Nahuatzen entre 0.15 y 0.24 ppm. 

La concentración de Pb ingerida en el organismo humano puede ocasionar severos 
dafios a la salud, por esta razón la OMS y múltiples normas internacionales y naciones 
establecen una cantidad máximo permisible. La norma mexicana NOM-250-SSA1-2014, 
establece 0.01 mg/L (ppm) como máximo permisible contenida en agua de consumo 
humano. Si se observan los resultados de la gráfica anterior, estos muestran estar decenas 
de ppm por arriba de la norma establecida, el cual muestra que debe ser un tema de 
puntual atención por autoridades y organismos. 


As (Arsénico) 


El arsénico fue otro de los elementos al que se prestó puntual atención para su 
estudio y determinación en las manantiales y pozos profundos, en la figura 6, se puede 
observar la concentración encontrada en cada uno de los sitios muestreados. En los 
sitios de muestreo de la localidad de Paracho se encontró un rango promedio de 0.018 y 
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0.02 ppm. Para los de la comunidad e Cherán oscila entre 0.06 y 0.1 ppm y en los de la 
localidad de Nahuatzen están en 0.1 ppm. Todos los sitios de muestreo revelan resultados 
por arriba de la norma permisible NOM-127-SSA1-1994, NOM-250-SSA1-2014 y la OMS 
el cual establecen una cantidad máxima de 0.01 mg/L ó 0.01 ppm. Si bien, no se han 
reportado resultados sobre el tema del arsénico y plomo en estos sitios y sus alrededores, 
cabe mencionar que la existencia de dichos metales y otros metaloides de interés pueden 
explicarse por varios factores presentes. La principal puede ser la presencia geológica, 
por ser una zona volcánica, por otro lado, investigaciones previas han reportado que 
el arsénico está presente en la atmosfera debido a fuentes naturales como reacciones 
ambientales, actividad biológica, emisiones volcânicas e incluso provenientes por un amplio 
rango de actividades antropogénicas (Campos et al., 2007; Alam et al., 2014) incluyendo la 
agricultura. Incluso el arsénico puede estar presente en el aire, el agua y el suelo a través 
de tormentas de polvo y agua de escorrentia, por lo que la contaminación de arsénico está 
muy extendida debido a su fácil dispersión (Gillispie et al., 2015 
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Figura 6. Concentración promedio en ppm (mg/L) de As en sitios de muestreo. 


CONCLUSIONES 


La calidad del agua para consumo humano es de vital importancia debido a su 
alta demanda y como medio de reacciones bioquímicas en el metabolismo celular. Por lo 
anterior se debe tener presente las normas vigentes nacionales e internacionales, para 
tomar en cuenta los parâmetros permisibles en la calidad de agua de mesa que ayuden a 
mejorar la calidad de vida. 
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En base a los resultados obtenidos en esta investigación se puede concluir que: 

-Se observó que las variables de dureza, SDT, pH y Temperatura se encuentran 
dentro de la norma establecida, NOM-127-SSA1-1994, NOM-250-SSA1-2014, la OMS e 
incluso por los establecido por la agencia de protección ambiental de los Estados Unidos 
de América. A diferencia de un sitio de muestreo en la comunidad de Cheran, el punto 
denominado manantial de Cumitzaro, la cantidad de dureza esta al rededor los 600 ppm y 
SDT está sobre los 400 ppm. 

-Por el contrario, para el caso del Pb y As, éstas se encuentran por arriba de la 
norma NOM-127-SSA1-1994, NOM-250-SSA1-2014. 
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RESUMEN: En este trabajo, se presenta el disefio 
de un molino compacto para triturar envases de 
polietileno tereftalato (PET), con el objetivo de 
que en una siguiente etapa se pueda implementar 
de manera física y hacer frente a la problemática 
que se tiene por la generación y uso excesivo de 
este tipo de envases, que si bien es cierto ofrece 
una gran ventaja para el envasado de bebidas, 
cuando no se les da un buen tratamiento se 
convierten en residuos que constituyen focos 
importantes de contaminación. Dentro del 
Instituto Tecnológico Superior P'urhépecha 
(ITSP), se observa la presencia de estos residuos 
en un promedio de 5Kg diarios, a pesar de la 
implementación de la normativa NOM-1S014000, 
por tal motivo se pensó en una solución para dar 
un tratamiento integral a estos residuos, siendo 
la primera etapa el disefo del molino, debido 
a que en la actualidad solo existen molinos o 
trituradoras de gran escala cuya adquisición es 
costosa. Por otro lado, en el mercado existen 
gran variedad de paquetes computacionales para 
el dibujo de piezas mecánicas, como por ejemplo 
el Autocad, sin embargo se optó por el uso del 
software SolidWorks, porque además de que 
permitió el dibujo de las piezas mecânicas de una 
manera relativamente sencilla, también brindó la 
opción de realizar un ensamble entre ellas, así 
como realización de cálculos de elemento finito, 
observando de esta manera las condiciones de 
operación del molino, permitiendo así contar 
con una validación del diseno, lo cual servirá 
para su implementación física. Los materiales 
utilizados en el disenfo, son de uso convencional 
y de fácil adquisición y por su tamafo puede 
resultar económicamente accesible para su 
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implementación dentro del ITSP o bien en alguna otra institución educativa pública o privada. 
PALABRAS CLAVE: Solidworks, Molino, Pet, Diseho, Simulación. 


ABSTRACT: In this work, the design of a compact mill to crush polyethylene terephthalate 
(PET) containers is presented, with the objective that in a next stage it can be implemented 
physically and face the problems that exist due to the generation and excessive use of this 
type of packaging, which, although it is true, offers a great advantage for the packaging of 
beverages, when they are not properly treated, they become waste that constitutes major 
sources of contamination. Within of Instituto Tecnológico Superior P'urhépecha (ITSP), 
the presence of these residues is observed in an average of 5Kg per day, despite the 
implementation of the NOM-IS014000 regulation, for this reason a solution was thought to 
give a treatment integral to these residues, the first stage being the design of the mill, since 
at present there are only large-scale mills or crushers whose acquisition is expensive. On the 
other hand, in the market there are a great variety of computer packages for the drawing of 
mechanical parts, such as Autocad, however, the use of SolidWorks software was chosen, 
because in addition to allowing the drawing of mechanical parts relatively simple way, it also 
offered the option of making an assembly between them, as well as carrying out finite element 
calculations, thus observing the operating conditions of the mill, thus allowing a design 
validation, which will serve for its implementation physical. The materials used in the design 
are of conventional use and easy to acquire, and due to their size, they can be economically 
accessible for implementation within the ITSP or in some other public or private educational 
Institution. 

KEYWORDS: SolidWorks, Mill, Pet, Design, Simulation. 


INTRODUCCIÓN 


La presente investigación se refiere al diseho de un molino compacio para triturar 
envases de polietileno tereftalato (PET). Entendiendo el diseão en ingenieria como 
“el proceso de aplicación de varias técnicas y principios científicos, con la finalidad de 
definir un dispositivo, un proceso o un sistema, con el detalle suficiente que permita su 
realización”[1]. La simulación en el software SolidWorks permite “someter fácilmente sus 
disefios a las mismas condiciones a las que se expondrán en el mundo real, permitiendo 
mejorar la calidad de sus productos, optimizar las geometrias y materiales, y además 
reducir costos tanto en tiempo como en el ahorro de prototipo para pruebas reales”[2]. El 
tereftalato de polietileno (PET) es un tipo de plástico muy usado en envases de bebidas y 
textiles [3]. Estos envases de PET son los que se presentan en la mayoria de las bebidas 
azucaradas que tienen una gran demanda en el mercado. El problema del plástico PET es 
que su producción es muy alta, y su reciclaje o reutilización se considera demasiada baja. 

En la actualidad sólo existen molinos industriales de gran capacidad y en 
consecuencia son de gran tamaho y costosos. Se pretende que el molino sea de bajo costo, 
comparado con los molinos comerciales, además de que los materiales utilizados serán 


accesibles y fácil de conseguir en el mercado, así como las partes de la máquina tienen una 
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forma geométrica sencilla y fácil de construir. 

El objetivo principal de este trabajo es la simulación del disefio mecánico de un molino 
de bajo costo para triturar envases o botellas de PET, provenientes del Instituto Tecnológico 
Superior P'urhépecha (ITSP). Debido a que, en la institución, se genera gran cantidad de 
botes de plástico que provienen de productos alimentícios (agua, jugos, refrescos, etc.) 
una vez que se han utilizado, los depositan en los botes de separación para su posterior 
venta o en su defecto tirarlos hacia a la basura. Al igual que en el instituto existen otras 
instituciones con el mismo problema anteriormente mencionado, el desperdício del PET. 
Para lograr el objetivo se recurren a las siguientes etapas del diseho; caracterización de 
las propiedades mecánicas del PET; propuestas de disehos; realizar el disefio más óptimo 
y crear la simulación en el software SolidWorks. Como el proyecto es ambicioso quedaria 


pendiente la fabricación de la máquina y la evaluación de ésta en otro momento. 


MATERIALES Y MÉTODOS 


La metodologia que se empleó para esta investigación fue el proceso de disefio de 
Robert L. Norton [1]. Este proceso consta de diez pasos que se muestran en la tabla 1. En 
este trabajo solo se desarrolló hasta el proceso 8. El proceso 9 y 10 se realizará en otra 
etapa de la investigación, que no está contemplada en el presente trabajo. 


EE ooo 
Tabla 1 El proceso de disefio 


Identificación de la necesidad 
Investigación preliminar 
Planteamiento de la meta 
Especificaciones de funcionamiento 
Sintesis 

Análisis 

Selección de la solución 


Disefio detallado 


Prototipos y pruebas 
Producción 


Õ 


Tabla 1. El proceso de disefo [1]. 


El paso inicial, identificación de la necesidad, normalmente consiste en un 
planteamiento del problema vago y escueto. Para este caso fue el no aprovechamiento de 
los envases de PET, provenientes de la institución. 


En el paso dos, investigación preliminar, se definió y se entendió por completo el 
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problema a resolver, así como también se investigó de los molinos actuales en el mercado, 
y se Ilegó a la conclusión de que solo existen de gran tamahio y en consecuencia son 
costosos. 

En el paso tres, planteamiento de la meta, se definió que la máquina debia de ser 
de medidas compactas, que los materiales se pudieran conseguir fácilmente en el mercado 
y como consecuencia se obtiene un molino cuyas características no se encuentran en los 
molinos del mercado actual. 

En el paso cuatro, especificaciones del funcionamiento, se define lo que el sistema 
debe hacer. Para esta investigación fue que la máquina entregue como producto final 
escamas de PET, para su posterior utilización. 

En el paso cinco, síntesis, es dónde se encontraron los enfoques de disefio 
alternativos que fueron posibles. Es donde se generó el mayor número de soluciones 
creativas. También se denomina el paso de ideación e invención — aqui se generó el mayor 
número viable de soluciones creativas. 

En el paso seis, análisis, se aplicaron técnicas de análisis por medio del elemento 
finito en el software SolidWorks, principalmente la Teoria de Von Mises-Hencky o de 
energia de distorsión, que es una magnitud física proporcional a la energia de distorsión. 
En ingenieria estructural se usa en el contexto de las teorias de fallo como indicador de un 
buen disefo para materiales dúctiles. 

La tensión de Von Mises puede calcularse a partir de las tensiones principales del 


tensor tensión en un punto de un sólido deformable, mediante la expresión: 


U = = [of + 0; + 0; — Z2v(g,0, + 0,203 + 0103)] (1) 


Donde: 

U, es la densidad de energia de deformación. 

E, es el módulo de Young o elasticidad. 

O, 0,, 0,, son los esfuerzos principales. 

ves la razón de Poisson. 

El esfuerzo efectivo de Von Mises o' se define como el esfuerzo de tensión uniaxial 
que crearíia la misma energia de distorsión que la combinación real de los esfuerzos 
aplicados [1]. 

o =V0) to +O — O; — 0,03 — 0,03 (2) 

La ecuación (2) también, se puede escribir como. 

S,=V0 +o;+o;-— 00, — 0,03 — 6,0; (3) 


a 


Donde: 
S, es el esfuerzo de fluencia del material. 
Para el esfuerzo bidimensional, 0,=0 y la ecuación (3) se reduce a 
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S, =V0Ê—- 00; +o; (4) 


” 


La ecuación bidimensional de la energia de distorsión (4) describe una elipse, la 
cual, al graficarse sobre los ejes 0,, 0,, queda como la de la figura 1. El interior de esta 
elipse define la región segura contra la fluencia bajo carga estática, para los esfuerzos 
biaxiales combinados [1]. Bajo esta teoria es que se analizaron las cuchillas o cortadores 


del molino, para no sufrir roptura. 






esfuerzo 5, 3 


” 
— Sy=0] —0/03 +03 =| 


se supone que Gs es 
el estucrzo igual a cero 


para torsión pura 
Sys= 05775, 


-1.5 -1.0 -05 0 05 10 1.5 esiuerzo / 5, 


Figura 1. Elipse normalizada de la energia de distorsión bidimensional, para la resistencia a la fluencia 
del material [1]. 


En el paso siete, selección de la solución, cuando el análisis técnico indica que se 
tienen algunos disefios potencialmente factibles, se debe seleccionar el mejor disponible 
para un diseho detallado, creación de prototipo y pruebas. El proceso de selección casi 
siempre implica un análisis comparativo de las soluciones de disefio disponibles [4]. Se 
realizó la comparación de diferentes disehos de cuchillas y se seleccionaron las más 
optimas en cuanto a disehfo y su resistencia a la falla. 

En el paso ocho, diseho detallado. Este paso en general incluye la creación de 
un conjunto completo de dibujos de ensamble detallados, o archivos de diseho asistido 
por computadora (CAD), por cada pieza utilizada en el diseho. Cada dibujo detallado 
debe especificar todas las dimensiones y las especificaciones de material necesarias para 
fabricar la pieza. Con estos dibujos (o archivos CAD) se necesita construir un modelo (o 
modelos) de prototipos para experimentos físicos [4]. Aquí se definieron las especificaciones 
técnicas del molino compacto, por ejemplo, la potencia del motor eléctrico, la transmisión de 
potencia por medio de poleas, disehio de poleas, tipo de bandas, el disefo de las cuchillas 
en el interior del molino, el tipo de material de nuestras cuchillas, y el tipo de base para 


sentar el molino. 
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DISCUSIÓN DE RESULTADOS 


Partes de la máquina del molino 


A continuación, se muestran las principales partes de la máquina del molino para 
triturar PET. La figura 2a, representa la vista frontal del juego de cuchillas, estas forman una 
figura tipo rehilete, el diâmetro exterior es de 123 mm. El cuadro en el interior es de lado de 
25 mm. El ângulo de la cuchilla forma 45º, debido a que, a ese valor, es cuando tenemos 
el cortante máximo. El espesor de las cuchillas es de 6.4 mm. Existe un separador entre 
cada cuchilla que tiene un diâmetro de 50 mm y en su interior tiene también un cuadrado, 
que mide cada lado de 25 mm. En la figura 2b, se muestra la vista superior donde van 
montadas las cuchillas y los separadores, en total tenemos 32 cuchillas y 33 separadores, 
como consecuencia tiene una longitud de 363.20 mm. Por último, la figura 2c, representa la 


vista isométrica del juego de cuchillas principales. 





Figura 2. Juego de cuchillas principales. 


La figura 3a, representa la vista frontal del juego de cuchillas base, tenemos una 
longitud de 148 mm y una altura de 50 mm. El espesor de las cuchillas es de 6.4 mm. 
Existe un separador entre cada cuchilla que tiene un espesor de 6.4 mm. En la figura 3b, se 
muestra la vista superior donde van montadas las cuchillas base y los separadores, en total 
tenemos 33 cuchillas y 32 separadores, como consecuencia tiene una longitud de 366.40 
mm. Por último, la figura 3c, representa la vista isométrica del juego de cuchillas base. 
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Figura 3. Juego de cuchillas base. 


La figura 4a, represenia la vista frontal de la caja cortadora, tenemos una dimensión 
de largo de 196 mm por 132 mm de alto. En la figura 4b, se muestra la vista superior de 
la caja junto con el eje de transmisión, este tiene una longitud de 789 mm. Por último, 
tenemos la figura 4c, que nos representa la vista isométrica de la caja cortadora y el eje de 


transmisión. 





Figura 4. Armado de la caja cortadora. 


En la figura 5, tenemos la base o mesa de soporte, dónde está montada la caja 
cortadora. La figura 5a, nos muestra la vista frontal, con dimensiones de largo 400 mm y 
de alto 1000 mm. En la figura 5b, se muestra la vista lateral de la mesa de soporte, con 
dimensiones de 850 mm por 1000 mm. Por último, tenemos la figura 5c, que nos representa 


la vista isométrica de la mesa de soporte. 
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Figura 5. Mesa de soporte. 


En la figura 6, tenemos las vistas principales representativas del molino compacto. 
La figura 6a, representa la vista posterior, donde se puede observar el conjunto de poleas, 
chumaceras, tolva de seguridad, mesa de soporte y las bandas de transmisión, para el 
funcionamiento de la máquina. En la figura 6b, se observa la vista superior y por último la 
figura 6c, muestra la vista isométrica del molino compacto para triturar PET. 


a) b) 


Figura 6. Molino compacto para triturar PET. 





Cc) 


Material 

El material que se utilizará para construir las cuchillas principales del molino 
triturador de PET, es el acero inoxidable 316L, debido a que se realizó un estudio de análisis 
electroquímico de diferentes aceros inoxidables 416, 304 y 316L, siendo éste el que tuvo 


mejor comportamiento. 
A continuación se muestra la comparación del comportamiento anódico de los 


aceros inoxidables 416, 304 y 316L. 
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Figura 7. Potencial V vs SCE. 


La figura 7, muestra las diferentes curvas del comportamiento anódico de tres 
aceros inoxidables, 416, 304 y 316L. Donde se observa el potencial aplicado en muestras 
de diferentes aceros sumergidos en una solución de H20-H2504, pH 1.0, Temperatura 
ambiente 20ºC. También es posible observar la densidad de corriente y las zonas activo- 
pasivo-transpasivo de las tres curvas. 


Generación del Código “G” 


Se realizó el código G de la cuchilla principal por medio del software Aspire. 
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Figura 8. Generación del código G, en el software Aspire. 


Ciencias exactas y de la tierra: Observación, formulación y predicción Capítulo 3 





Este es el código G de la cuchilla principal. 


NOG90 

N10(End Mill (0.25 inch)) 

N20M6T1 

N30M3 S12000 

N40GO X-0.3975 Y0.5520 Z0.2362 

N50G1 Z-0.0984 F30.0 

N60G3 X-0.0490 Y-0.3578 10.9642 J0.5520 
F100.0 

N70GO Z0.2362 

N80GO X-0.3975 Y0.5520 

N90G1 Z-0.1969 F30.0 

N100G3 X-0.0490 Y-0.3578 10.9642 J0.5520 
F100.0 

N110GO Z0.2362 

N120GO X-0.7230 Y-1.2159 

N130G1 Z-0.0984 F30.0 

N140G1 X-0.4057 Y-1.4933 F100.0 
N150GO Z0.2362 

N160GO X-0.7230 Y-1.2159 

N170G1 Z-0.1969 F30.0 

N180G1 X-0.4057 Y-1.4933 F100.0 
N190GO Z0.2362 

N200GO X-0.3853 Y-1.5052 

N210G1 Z-0.0984 F30.0 

N220G3 X2.3258 Y0.5562 10.2910 J0.4189 
F100.0 

N230GO Z0.2362 

N240GO0 X-0.3853 Y-1.5052 

N250G1 Z-0.1969 F30.0 

N260G3 X2.3258 Y0.5562 10.2910 J0.4189 
F100.0 


N270GO Z0.2362 

N280GO0 X2.3259 Y0.5520 

N290G1 Z-0.0984 F30.0 

N300G3 X1.9775 Y1.4617 10.9642 JO.5520 
F100.0 

N310GO Z0.2362 

N320G0 X2.3259 Y0.5520 

N330G1 Z-0.1969 F30.0 

N340G3 X1.9775 Y1.4617 10.9642 JO.5520 
F100.0 

N350GO0 Z0.2362 

N360GO0 X2.6514 Y2.3198 

N370G1 Z-0.0984 F30.0 

N380G1 X2.3341 Y2.5972 F100.0 

N390G0 Z0.2362 

N400GO X2.6514 Y2.3198 

N410G1 Z-0.1969 F30.0 

N420G1 X2.3341 Y2.5972 F100.0 

N430G0 Z0.2362 

N440GO X2.3137 Y2.6091 

N450G1 Z-0.0984 F30.0 

N460G3 X-0.3974 Y0.5477 11.6375 JO.6850 
F100.0 

N470GO Z0.2362 

N480GO0 X2.3137 Y2.6091 

N490G1 Z-0.1969 F30.0 

N500G3 X-0.3974 Y0.5477 11.6375 JO.6850 
F100.0 

N510GO Z0.2362 

N520M30 





Tabla 2. Código G de la cuchilla principal. 


A continuación, se muestran los análisis que se realizaron por medio del elemento 
finito, mediante el software SolidWorks. En la figura 9. Se encuentra el análisis de Tensión 
de Von Mises. Se puede observar que las cuchillas no están en riesgo de fracturarse, 
debido a que el máximo esfuerzo que tienen es de 1.586E8 (N/m?) y el límite elástico del 
material de la cuchilla, que es de acero inoxidable 316L es de 1.7E8(N/m?). Se puede decir 
que se tiene un factor de seguridad de 6.7%, usando las cuchillas en su condición. 
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Mombre del modedo:cuchiliaotrad 
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Figura 9. Análisis de Tensión de Von Mises. 


En la figura 10. Se muestra el análisis de desplazamientos resultantes, se puede 
observar que el máximo desplazamiento resultante, se encuentra en la parte del filo de 
la cuchilla, que es de 0.0391547 mm, pero este desplazamiento no es significativo para 
nuestro diseho, debido a que no alcanza ni un milímetro, por lo tanto no afectaría con las 
tolerancias establecidas en el disefio. 
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Figura 10. Análisis de desplazamientos resultantes. 
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CONCLUSIONES 


El uso de materiales poliméricos en la sociedad actual a crecido en forma exponencial, 
debido a la producción y uso. Sin embargo la utilización de polímeros, ha Ilevado a la 
sociedad a generar un gran problema de contaminación, particularmente el Polietileno de 
Teraftalato (PET). Para ello se vio la necesidad de disefiar y modelar un molino compacto 
con el programa SolidWorks. 

Se desarrollaron todos los elementos de máquina del molino a través del programa 
SolidWorks, desde el inicio hasta el ensamble de cada uno de ellos. Se puede observar 
que no hay interferencia en cuanto a las partes que componen la máquina, lo que significa 
que, en la construcción en físico, no deberia tener ningún problema en su funcionamiento. 

Se obtuvo una gráfica del comportamiento anódico de diferentes aceros inoxidables, 
para la fabricación de las cuchillas, se determinó que, el acero 316L es el mejor material, 
debido a su resistencia electroquímica. 

Se generó el código “G” para la fabricación de las cuchillas, en el software Aspire. 
Una vez que se realicen en físico, las cuales se van a realizar en una fresadora de control 
numérico. Los código se adaptaron para trabajar en una fresadora ExperMill VMC-0600. 

Se realizó el análisis del elemento finito de las cuchillas en el software SolidWork, 
los análisis que se realizaron fueron el de Tensión de Von Mises y de Desplazamientos 
resultantes, en ninguno de los dos análisis existieron problemas, debido a que no se superó 
el límite elástico del material y también el desplazamiento resultante fue menos de medio 
milímetro. 

En la actualidad podemos realizar el diseno y la simulación de una máquina, antes 
de realizarla en físico, esto nos ahorra bastante tiempo y dinero. Se puede asegurar que 
una vez fabricada, funcionará adecuadamente debido a que ya se realizó la simulación con 


los modelos o teorías adecuadas. 
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RESUMEN: Se describe de forma breve la 
interacción entre ondas electromagnéticas, 
mecânicas y sonido debido a una perturbación 
de un cuerpo por transmisión de energia. 
Las funciones para diferentes condiciones 
del experimento se muestran en gráficas. 
Solo algunos datos aparecen en las tablas. 
Se obtuvieron cientos de datos durante el 
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experimento, debido a los tiempos para cada 
prueba, con un máximo de un minuto entre 
cada prueba de 30 segundos, para no agotar 
la elasticidad de la membrana del micrófono 
receptor. Con la inclusión de la ley de Hooke para 
una membrana delgada que vibra sin deformarse, 
se llegó a una modificación de la ecuación 
diferencial del modelo de cuerno acústico cuya 
solución y desarrollo estân comentados en 
algunas secciones. La perspectiva de la utilidad 
de la fotoacústica unida a su porqué, su cómo y 
aplicación, se muestran en este trabajo a través 
de correlacionar la simulación, el experimento y 
el análisis matemático de un modelo en particular. 
PALABRAS CLAVE: Efecto fotoacústico; hoja 
metálica; densidad; grosor; Ley de Hook; láser 
rojo. 


ABSTRACT: This work shows math analysis of 
interaction for electromagnetic, mechanics and 
sound waves due to an energy transmission 
on a thin plate like body. Graphs correspond to 
some functions about different experimental 
conditions. Tables present some of hundreds 
of measured data. Those values come from 
multiple experiments of 30 seconds each, with a 
relaxing time of one minute. More than that use 
generate a loss of homogeneous elastic behave 
on microphone. Based on varying-cross sectional 
area rod differential equation, with known answer, 
emerges a new differential equation including a 
Hooke's law type element. According to use of 
thin metal plate sheet. This paper is a general 
perspective about photoacoustic effect use on 
industry, how to improve it and advantages. 
All through simulation of differential equation, 
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experimental data and math analysis for our particular model. 
KEYWORDS: Photoacoustic effect; metal plate sheet; density; thickness; Hooke Law; red 
light laser. 


11 INTRODUCCIÓN 


Las perturbaciones y ondas se encuentran en los fenómenos de la naturaleza. Se 
hace necesario conocer la forma en que estas interaciúan. La fotoacústica está generando 
un mayor interés debido a su utilidad para el análisis de propiedades mecano-elásticas, 
mecánicas y visco-elásticas de materiales y en comunicaciones. 

Se tiene como tendencia el uso de láseres de alta potencia (> 1W) y láseres 
concentrados, como se reporta en la literatura. En contraste con esto, se trabajó con un 
diodo láser de baja potencia (50 mW) para generar el fenómeno fotoacústico. Durante las 
pruebas se utilizó un sistema dual de micrófonos en donde uno genera una sefial de 40 kHz 
como referencia y el otro, mide el efecto producido en el objeto por interacción del láser. 
Diodo láser de una longitud de onda de 650 nm. Situado por encima de los micrófonos, al 
centro y apuntando hacia la placa metálica con un ângulo no mayor a 15 grados. 

Se aprovechó la simetria en cuanto al arreglo. Considerando que la luz genera una 
ligera presión sobre la placa [LANDAU Y LIFSHITZ, 1963, RANDALL, 2005, ELMORE, 
HEALD, 1985], una compresión que después se interrumpía. Este ciclo de presión-relajación 
descrito por ley de Hooke genera una onda sónica de baja intensidad que provoca una 


alteración en la referencia (fase y amplitud). 


2 | ANÁLISIS MATEMÁTICO 


Usando sistema dual de micrófonos (forma cilíndrica, modelo HC-SRO4 
[HARDWARELIBRE]) con el rayo sónico con la sefial de referencia ultrasónica de 40 kHz. 
Para el análisis de la respuesta del efecto fotoacústico se utiliza un sistema coordinado 
cuyo origen está en el centro simétrico del HC-SR04. Ambos micrófonos tienen un disefio 
con inclinación de 15 grados orientados hacia el eje central. Equidistantes con el punto de 


origen. 
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+ air (medium) 


PAR: 
sensor control f a» 








T = rpom temperature |) 
Figura 1. Bosquejo para determinar el diagrama de fuerzas y la dinámica de la onda. 


Para la onda viajera de ultrasonido, se considera que las variaciones de densidad 
son suficientemente pequefas en el modelo matemático, permanece constante. Se hace 
el análisis usando de referencia el eje principal. Para las perturbaciones en el medio, 
hay pequefios corrimientos de fase y cambios de intensidad. El HC-SRO4 permanece en 
funcionamiento durante las pruebas, de manera que se desprecia el decaimiento de la 
sefial acústica en el tiempo. 


Las consideraciones pertinentes para el análisis del fenómeno son: 
* | No hay gradiente de presión en la zona de las pruebas. 


* | No hay variaciones de temperatura debido al fenómeno [WILEY, RANDALL, 
2005, NAVARRO). 


* | Elgrosor de placa es despreciable con respecto a su extensión: se considera al 
cuerpo como una membrana vibrante. 


. Se consideran las densidades del cuerpo como del medio, constantes. 
* - Setrabaja bajo el supuesto de un grosor promedio de la membrana. 


Con el diagrama de fuerza se genera una ecuación a partir de la tercera ley de 
Newton, con la acción de la perturbación del aire, la respuesta de la placa delgada y la 
resistencia de micro gotas de agua en el medio. Se considera que el aire y vapor de agua 


funcionan como fuerzas generadoras de la perturbación: 


p=vysl, as e E 
“rr cdi dr (1) 

Cuyos parâmetros son: módulo de Young (Y), sección transversal (9), viscosidad 
dinâmica del líquido (u, agua para este caso), frecuencia angular de vibración respuesta 


(ww), la función que describe la perturbación (1). 
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Condición de equilibrio: Se considera que el aire reacciona a las compresiones 
por la membrana. Esta fuerza de reacción se relaciona con una pseudo aceleración, y un 


término que describe el efecto de ley de Hooke: 


dn | | 
F = PoST 572 — (—Y dr cos wt) (2) 


No se trata de una interacción cuerpo-a-cuerpo per se. En cambio, se analiza el 
cambio de las fuerzas en la zona de contacto, correspondiente a la sección transversal de 


la perturbación acústica. En un símil a la tercera Ley de Newton: 


Si yvsa prt ara O e yidesoi 3 
gra AS qo = os gpa O soênio 1) 


Dicha ecuación corresponde a los cambios tangenciales de la perturbación. Con 
algebra, despeje y renombrando constantes, se Ilega a la ecuación diferencial: 
dn 1 9“n BD 


dr? “(Lra)iot? '* (L+a) 





Cos wt (4) 


Donde las constantes y parâmetros de utilidad son: 


-  c= VY/po velocidad de la perturbación, descrita por el módulo de Young y 
densidad del aire a 20 ºC 


: a = w/a factor de atenuación de la viscosidad 

“ a=Yu 

*- B=Y/Yconstante sin unidades, establece el límite útil del modelo 
- D=da/S 


Se propone una solución a través del método de separación de variables: 





n(r,t) = fO)T(L) = f(r) cos wt (5) 
Y con su sustitución en la ecuación diferencial se tiene: 

d?f ww? Pro. BD 

—— —— 2 —— h E = ————— 

de ara “qe O ra) (6) 


Considerando el término & = (Y'/Y), como un límite a dicha teoria establecida en 
el trabajo solo se puede tener una variación de unos cuantos ordenes de magnitud, entre 
un módulo de Young y el otro. Si Y>>Y, B tiende a infinito, por consiguiente, también la 
función, algo fisicamente imposible. Si Y>>Y, £ tiende a cero, entonces el pulso acústico 
solo podria existir si la amplitud inicial es muy grande, es decir, para una fuente láser de 
alta potencia. 

Un análisis de unidades del término de ley de Hooke, F, = -Y“dr, muestra que Y'd 


tiene unidades de constante de resorte. La expresión de amortiguamiento, (1+a.) [MILNE- 
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THOMSON, 2011, MORFEY, SOROKIN, GABARD, 2012], no se puede despreciar incluso 
si el nivel de humedad relativa es pequefo. Se busca demostrar que el análisis matemático 
es válido para condiciones con un nivel de humedad diverso. 
La ecuación diferencial homogénea de una onda en una barra de densidad constante 
y sección transversal variable [ELMORE, HEALRD, 1985], muestra una similitud con la 
establecida previamente: 
dºf dºf k* 
-—+hif=>+———f =0 7 
dr* 1f dr* (1 To)! Ea 
El análisis físico requiere la solución senoidal. Para el caso de la ecuación diferencial 
no homogénea, se recurre a la integral de Green [BUTKOV, 1973]: 
f0)= | Gr;Drar= | sink(£-nag=- [sinkr- ds (8) 
Fo Tn [] 


Cuya solución (para el caso en que r =0) no tiene componente complejo: 








“ABD 
fo=DDa-cstr) O) 
p= NE (10) 


Se requiere una corrección para coincidir con los datos experimentales. Por 
condiciones del sistema de toma de datos, la amplitud inicial se considera constante, con 
A,=9V. Se tiene que la función solución consta de dos partes, una para el caso en que no 


hay oscilación armónica en la membrana, como se muestra: 
n(r,t) = Ag(1l — coskyr) cos wt (11) 


Un cuerpo sólido perfecto, se conoce como cuerpo rígido. Esta placa, a diferencia, 
reacciona con una vibración en la zona cercana al impacto de la sehal de referencia. Una 
vez que el HC-SRO4 manda otro pulso, se tiene un tiempo breve para que la membrana 
se relaje antes de volverse a compactar. Estas vibraciones armónicas dependen de las 
características del medio y la placa. 

La resonancia se presenta en algunos materiales cuando su extensión es 
debidamente larga o cuando su grosor es suficientemente pequefio. En este modelo es 
necesario considerar dicha resonancia. Se hicieron pruebas con placas metálicas delgadas 
(orden de mm). Considerando dicha resonancia cual un repetidor de frecuencia angular 
como parte del vector de desplazamiento. Los valores posibles para este factor son enteros 
positivos. La resonancia se puede detectar como resultado del camino de regreso, o parte 
del recorrido completo. 

Supóngase que el factor de resonancia es: 
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LES E a 145 TO dq 
«avl+a da cyvl ta 
Con lo que se obtiene una función más precisa: 


— n2k 
=A|1+ l-—- cosk,r)coswt (14 
p= (14 E Eta ent dio 





Este factor de la resonancia (o repetidor de la amplitud) no desaparece si se sustituye 
la nueva función solución en la ecuación diferencial no homogénea. Donde dicho factor, 
predice una variabilidad de la amplitud medida. 

Una vez establecida la solución, es necesario establecer la relación con las 
propiedades a medir. Si el cuerpo en estudio fuese un cuerpo rígido súper masivo, la sefial 
de referencia haría el recorrido sin tener un cambio en amplitud. La placa en cuestión es 
un cuerpo que reacciona a la sehial incidente generando una pequefia vibración en la zona, 
donde esta cambia la uniformidad de dicha sefial recibida. Cambia la amplitud (aumentando 
o disminuyendo) de acuerdo a la vibración de la membrana. 

Se puede aproximar que la distancia recorrida por la sefal es equivalente a la 
separación del HC-SRO4 con la placa y el grosor de la placa. Esta variación se da en 
términos de la sefial original, A, y la medida, A,; 


À cvl+a s (e — +=) (15) 
1 =——— cos —— 
Zu) Ao sia 


En el tiempo cero, el sistema permanece en equilibrio. Una vez que se enciende el 
HC-SRO4, se genera una sefial ante la cual, la membrana intenta oponerse en los intervalos 
de tiempo que no hay sefial enviada (1 s entre pulsos). Se puede ver como dos resortes 
interactuando entre ellos (el medio y la placa delgada); 


FSAr 





Hr=d cos wt = Faire = Feuerpo = (Y'd)d cos wt (16) 


Si consideramos que Ar=1 unidad, gracias a que la sefial permanece casi constante, 
se reescribe la ecuación (16) de forma que se pueda despejar el módulo de Young 


correspondiente a la placa delgada; 


, Nr=df'S 
V = FE (17) 
Para poder medir la densidad de la placa, se recurre primero a su expresión 
simplificada: 
EM na = [dl 18 


Desde el comienzo del desarrollo se estableció la interacción de dinâmica de fuerzas 
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con apoyo de las leyes de Newton. El término correspondiente a la aceleración viene de la 

idea que se trata de una barra que vibra (cilindro de aire) [7, 8]. Relacionando este concepto 

con la ecuación (18) se puede medir la densidad del cuerpo: 
m Ffja. (Y' dir cos wt 


Pg” as fas 
IE 


(19) 





Donde W=A, cos krcos wt, la densidad se despeja como: 


di 


“ SwlAgcosk;r (20) 


Pp = 


Si la sehal medida por el HC-SRO4 después del recorrido, tiene una relación con 
la frecuencia natural, se puede decir que el corrimiento de fase está relacionado con esta 
senial. Se puede conocer la frecuencia de resonancia natural, en el caso que se conozcan 


las características del cuerpo: 





(d) = À (1 Pl id = A (21) 
fld) = Ag ic, pe 6 =) - 2 el 
x cvl+ta o. (Ag— Às 

d=w- 7 es (DS) (22) 


O 


31 ENFOQUE NUMÉRICO 


Se usó el método de resolución de ecuaciones numérica de Runge-Kutta de 
cuarto orden, para la ecuación diferencial radial, suponiendo una función temporal del 
tipo trigonométrica como se expuso anteriormente. La finalidad de conocer la solución 
radial pura es poder determinar cuan rápido decae la onda generada para las mediciones 
experimentales. 

Para el código en MatLab (Licencia 40673431), se establece como aclaración lo 


siguiente: 


À k = 


tal 


cvl+o 





vector de desplazamiento, relacionado con la densidad. 


Para el primero conjunto de simulaciones, se busca conocer que tan pronunciado 
es el cambio de amplitud de la onda generada en la placa, de acuerdo a la ecuación (6). 
Las condiciones pertinentes para la simulación y los parâmetros fijos se mencionan a 


continuación: 
* | Rango deradio:0 a 5 
* | M. de Young del aire [WIKIPEDIA]: 0.2 MPa 
*- M. de Young del acero [WIKIPEDIA]: 210 000 MPa 
* | Velocidad de sonido en el medio: 343 m/s (a 20 ºC) 


. Velocidad de onda inicial: O m/s 
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* | Grosor de placa: 0.0001, 1.5 
. Factor de reducción (alfa): 0.0001, 1, 15 
* | Frecuencia de pulso: 40 000 Hz 


Las soluciones para la ecuación diferencial no homogénea son exponenciales, 
trigonométricas y/o superposiciones de estas. Pueden ser obtenidas por diferentes métodos 
de resolución analítica o numérica. También, dependiendo de condiciones frontera o de 


INÍCIO. 


Onda acústica (medida) 








0 0.5 1 1.5 2 2.5 3 ER) 4 4.5 5 
r- distancia adimensional 


Figura 2. Función radial esperada para un grosor de 0.0001 


Onda acústica (medida) 
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Figura 3. Función radial esperada para un grosor de 1.5 


En las simulaciones realizadas, se obtuvieron funciones de tipo exponencial. Lo 
pronunciado en las gráficas cambia de acuerdo a los parâmetros de grosor y el valor de alfa 
(relacionado con el amortiguamiento debido de la humedad relativa). 


La sehal recibida en el segundo micrófono, es una nueva función de tiempo, con un 
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posible corrimiento de fase. La siguiente gráfica muestra el comportamiento de la sefial en 
términos del tiempo y posibles fases, para valores de simulación: alfa = 1, grosor = 0.0001. 


Onda acústica (medida) 


- onda 


ftr) 





0 0.5 1 Te 2 20 3 3.5 4 4.5 3 
r- distancia adimensional 


Fig. 5. Perturbación con su amplitud máxima para diferentes frecuencias de respuesta. 


La gráfica obtenida muestra una función que parece incrementar hasta el infinito, 
esto se debe a que la ecuación diferencial como la solución, son válidos para la zona del 
fenómeno fotoacústico y no se tienen condiciones de amortiguamiento debidas a la pérdida 
de homogeneidad de acuerdo al aumento de temperatura o de la fricción que el gradiente 
de temperatura usada y las condiciones del medio. En condiciones controladas, el modelo 
es una buena aproximación para sistema de control de calidad en una empresa para un 


medidor ubicado cerca de la línea de producción. 


41 RESULTADOS Y TRATAMIENTO DE DATOS 


Sin un circuito con filtrado de ruido apropiado, el sonido de fondo está presente en 
la sefal medida. Para cada intervalo de tiempo, la amplitud registrada tiene un ruido debido 
al diseho del circuito y ruido dispersivo debido a que la sehal ultrasónica se debilita a causa 
de la potencia. Aun si se considera que el gradiente de temperatura no forma parte del 
fenómeno simplificado estudiado. 

Para hacer una clasificación apropiada de los datos medidos, se debe tener 
primero un registro del ruido electrónico como el ruido de fondo del medio gaseoso. Sin 
la presencia de la fuente láser. Cuidando que cada registro de ruido como las mediciones 
experimentales, tuvieran la misma temperatura ambiente (295.15 K) y presión atmosférica 
(alrededor de 101 355 Pa). 

En el caso del registro de ruido como el registro de datos experimentales se 
calcularon los promedios temporales. A considerar en el fenómeno estudiado tres tipos de 
perturbaciones: el ultrasonido incidente en el cuerpo, la resonancia del cuerpo y el efecto 
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de ley de Hooke debido al láser. Posteriormente, con dichos registros se calculan los datos 
de los parâmetros medidos. 

De la ecuación (12) se obtienen las relaciones entre las amplitudes medidas como 
las de referencia, para el análisis de los datos. Dentro de las pruebas realizadas, se trabajó 
brevemente con cartón delgado debido a su importancia en la industria papelera y su 
aplicación para una empresa mexicana. Sin embargo, este último material es desconocido 
en sus componentes básicos. 

En la tabla siguiente, aparece el término “factor de corrección”, que es el coeficiente 
numérico requerido para los valores de amplitud en que se satura el sistema del micrófono. 
Del factor de repetición en la ecuación (15) se puede hacer un despeje de la sefial alterada 
con respecto a la superposición de ondas de la sehal base con el ruido de fondo y el 
ruido del HC-SRO4. Ahora, en la siguiente tabla se muestra el “% de interacción”, o bien, 
el nivel de resonancia (12) encontrado para una sefial medida y para ciertas placas. Y 
posteriormente se muestran los datos obtenidos experimentalmente con su error relativo, 


calculado como sigue: 


” | Preórico — 


experimental | 
Error = ao 


100% (23) 


Vreórico 


Amplitud Densidad ; 

Grosor real Factor de Amplitud sin e Módulo de Young 

(mm) dis ed li s corrección resonancia OE RO: reportado (MPa 
2 


mº) 


) 
7 850 
9.709551 [WIKIPEDIA] 210 000 [BUTKOV] 


0.650 2 06896552 2 06896552 7 850 210 000 
0.610 269230769 | 0732272 | 1.96003569 7 850 210 000 
0.406 10.222222 | 10.98408157 | 0.76185957 7 850 210 000 
0.305 6.78571429 | 7,3922721 0.53700671 7 850 210 000 


Tabla 1. Valores de referencia para las comparaciones y cálculo de error. 





Frecuencia angular 251 327.4123 Hz Velocidad del sonido 343 m/s 
viscosidad, u 1002 uPa*s (agua) 
Factor (vector de Factor (frecuencia para la Dciciscs Cain 
desplazamiento) resonancia de la onda) 
0.732272 502 654.8246 (aprox.) O 
2.929088 753 982.2369 
75 


1 005 309.649 66.66667 


| Frecuenciaangular — 

| viscosidadu 
n,m 

e 
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35.88133 2 010 619.298 85.71429 


o 9 59.31404 2 513 274.123 g8.88889 


Tabla 2. Algunos de los factores de resonancia que surgieron en los datos experimentales. 





Se muestra una breve descripción de cada factor involucrado en el análisis: 


* | Dejando de lado la composición atómica y estructura cristalina de los cuerpos, 
en este caso, de acuerdo a las dimensiones, se considera a la placa delgada 
un cuerpo unitario. Esta percibe la perturbación como una onda que viaja por 
dentro de ida y vuelta. Válido solo para ciertos rangos de tiempo de respuesta 
y de interacción, propios para cada material. Con esto, se obtiene el factor de 
vector de desplazamiento mencionado, como la cantidad de veces que la onda 
“rebota” dentro de la placa: 


á Z z 
n2ky = = = EL (24) 


vi+a cvl+a 








* | Enotros casos, si consideramos la estructura molecular y/o atómica, de forma 
que se puede tener una resonancia en la membrana. Cada partícula tiene su 
valor propio de vibración, y en la interacción de los átomos conformando el ob- 
jeto se tiene una frecuencia de resonancia propia para el material. Sin conocer 
los valores de resonancia del material, se puede suponer que la serial del HC- 
-SRO4 tiene un tiempo efectivo en la placa, para que se sume la serial original y 
se puede conocer este factor a través de: nw=(1+) w=w+lw. 


à Erosar Módulo de 
(sin % error Young % error densidad, p % error 
medido 
resonancia) (- 10" N/m?) 


1.900 | 9709551 | 1911222209 | 0.5906 | 0.4281949543 RO 1559.314219 ER 
1.800 | 3.580466 | 1.752098776 | 26612 0.2049461403 746.3316148 


0.650 | 2.06896552 | 0.644455345 | 0.8530 2.3/9856174 13.3265 8666.481344 10.4010 
0.610 | 1.96003569 | 0.625921228 | 2.6100 2.2094558132 7.30991 8210.196147 4.5885 





Tabla 3. Medición de los parámetros y cálculo del error relativo. RF = Fuera de Rango 
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Grosor real (mm) Amplitud (medida) Factor de Frecuencia de Factor de repetición 
p a E medida ERR 
0.406 10.222222 | 6590449 590449 - | -2985528171 552.8171 1. | AT5sB783D 
0.305 6.78571429 2.929088 - 502 667.3546 2.000049855 


0.203 2 894736842 - 709 162.0716 2 821666229 
0.152 20.27586207 18.306800 - 786 830.7529 3.130700059 
0.102 3.344262295 0.732272 - 1553 707.389 6.182005277 


Tabla 4. Análisis de datos para factores de repetición relacionados a valores de vibración. 





Con estos datos se puede concluir la solución para el análisis de propiedades: 





An(l—- cosk;r)coswt; | sinresonancia 
oo 
“= nºk mu 
n(r,t) Ay (1 + Ea (1 — Cos ===) coswl;  resonancia (2a) 
| a 4 l+a l+a 


5 | TRABAJOS FUTUROS 


El efecto de la ley de Hooke genera una perturbación que viaja por el cuerpo en todas 
las direcciones, dinâmica que fue omitida en el análisis. Dichas consideraciones se deben 
incluir para trabajos futuros, junto como la posible deformación del cuerpo o compresiones 
que requieran un tiempo grande para recuperar su forma. A futuro se planea usar el tensor 
de estrés y deformaciones en alguna forma simplificada que involucre cuerpos isotrópicos, 
pero con medios anisótropos. Los elementos matemáticos que no se han considerado 
sobre deformaciones que evolucionan con el tiempo o variaciones en sistemas casi reales 


de la densidad, son parte del siguiente escalón. 





Figura 5. Corte transverso de la estructura de la placa metálica. Representación básica del efecto de la 
ley de Hooke propagándose en las direcciones posibles. 


é En qué nivel se induce la resonancia en el efecto fotoacústico? Con la luz 
interactuando con el metal, ;toma importancia el efecto de skin depth? é Qué ocurre cuando 
hay una polarización del medio? Todas estas son cuestiones que se deben responder en 
trabajos futuros. Incluyendo el uso de un modelo matemático más completo que pueda 
prever ciertas condiciones como la falta de simetria, la resonancia o materiales que no son 


tan rígidos. 
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61 DISCUSIÓN 


La física permite conocer ciertas condiciones del efecto fotoacústico que no aparecen 
en estas simulaciones. La experiencia muestra que una sehal acústica eternamente 
creciente no existe. Conforme una onda acústica evoluciona en el tiempo, se genera 
un gradiente de la densidad y un cambio de la temperatura, disminuyendo la primera, 
dificultando la posibilidad de propagación de un sonido (incluso ultrasonido o infrasonido). 

No hay forma de estudiar el efecto fotoacústico en su totalidad con unas cuantas 
simulaciones sobre una ecuación diferencial homogénea. Por su complejidad requiere 
simulaciones y estudios de los efectos internos, el gradiente de temperatura de la placa, 
adherencia de la humedad, la importancia de los eigenvalores y las posibles deformaciones. 

Otros de las consideraciones necesarias, es el estudio de los tipos de ruido que 
están involucrados y como se puede analizar estos. Incluso ver si los tipos de ruido pueden 
ser eliminados sin perder información. Incluso considerar el ruido de Johnson como parte 
del fenómeno y no solo hacer pruebas con variaciones pequefas de temperatura. 

La evolución de la función en el tiempo y la distancia, nos da información sobre 
el sistema como algo libre de restricciones. Es de gran importancia conocer el valor de 
distancia límite del cambio significativo de la presión, alterando la densidad, la viscosidad 
o ambas. Con eso se puede tener un análisis no limitado a casos adiabático e isotérmico. 

Se demostró que los factores de resonancia y de repetición son importantes en 
la caracterización de metales por efecto fotoacústico. Datos experimentales muestran un 
error relativo considerablemente pequeho. Para conocer los valores precisos de n, m se 
requiere establecer un nuevo modelo matemático o un conjunto de experimentos en que se 
conozcan las propiedades de las placas. Para placas muy delgadas, el factor de resonancia 
es mayor. 
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RESUMEN: En esta investigación científica se 
estudió numéricamente el modelo matemático de 
Harderberg — Lefever. Este modelo es un sistema 
de ecuaciones diferenciales parciales no- 
lineales, el cual, reproduce una amplia gama de 
patrones observados con regiones limitadas en 
agua, incluida fase de puntos entremezclada con 
fase de desertificación con laberintos. Predice 
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las transiciones de suelo desnudo con poca 
precipitación a la vegetación homogénea con alta 
precipitación, pasando por estados intermedios 
de patrones de puntos, tripas y huecos. 
Utilizando estas predicciones, proponemos una 
explicación novedosa de los fenómenos de 
desertificación y un nuevo enfoque para clasificar 
la distribución de biomasa en forma de patrones 
espacio-temporales. Estos resultados muestran 
que, en escenarios climáticos futuros, el destino 
del sistema dependerá de las características 
y las adaptaciones de las plantas, así como, 
de la influencia que tienen las condiciones de 
humedad sobre ellas. 

PALABRAS CLAVE: Patrones de vegetación, 
desertificación, ecohidrologia. 


MODELING VEGETATION PATTERNS 
WITH NONLINEAR EQUATIONS USING 
THE HARDENBERG AND LEFEVER 
MODELS 


ABSTRACT: This paper proposes a new model 
to describe the vegetation distribution of biomass- 
desertification. This model is a system of nonlinear 
partial differential equations, which reproduces 
a wide range of patterns observed in water- 
limited regions, including points and intermixed 
phases of desertification with mazes. It predicts 
transitions from bare soil with low precipitation to 
homogeneous vegetation with high precipitation 
through intermediate states of point, gut, and gap 
patterns. Using these predictions, we propose a 
novel explanation of desertification phenomena 
and a new approach to classify the biomass by 
means of the morphology of their pattern. The 
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results obtained by the model show that, in future climate scenarios, the fate of the system 
will depend on the characteristics and improvements of the plants, as well as the influence of 
moisture conditions on them. 

KEYWORDS: Vegetation-patterns, desertification, ecohydrology. 


11 INTRODUCCIÓN 


El cambio climático es un tema que, durante aos, ha capturado la atención de todo ser 
pensante en el globo terráqueo. Los temas ambientales han sido ignorados o considerados 
secundarios, debido a que no se le ha dado una justa dimensión a la importancia de un medio 
ambiente sano; o bien, ingenuamente como especie humana conservamos la esperanza de 
que al solucionar los problemas sociales y económicos los ecosistemas naturales tendrían 
tiempo y capacidad de recuperarse por sí solos. Sin embargo, se enfrenta el hecho de que 
el deterioro ambiental amenaza seriamente el desarrollo actual y el futuro de diferentes 
ecosistemas [13]. La degradación del medio ambiente con el consecuente cambio climático 
es una bomba de tiempo que debe ser desactivada, si el ser humano no quiere desaparecer 
como especie de nuestro planeta. 

Se han sugerido una gran cantidad de hipótesis sobre el origen de formación de 
patrones de vegetación. Existe una gran variedad de modelos matemáticos que describen la 
formación de estos mismos o especificamente: la distribución vegetal de biomasa mediante 
las leyes de la Física y su comportamiento no lineal para modelar la interacción de las 
plantas con su medio ambiente. Este nuevo esfuerzo científico ha fructificado en nuevas 
áreas de la ciencia como la Ecohidrologia [1, 7] y otras herramientas como la Alometria 
[10]. Tal es el caso del modelo de Hardenberg [11], el cual predice las transiciones de 
suelo desnudo con poca precipitación y la vegetación homogénea con alta precipitación 
donde coexisten diferentes estados estables. O bien, el modelo de Lefever [17] donde los 
resultados modifican la interpretación común que asocia los patrones de vegetación y el 
proceso de desertificación local bajo inestabilidad de Turing, dando así una comprensión 
unificada de la autoorganización de la vegetación dentro del amplio contexto de las 
transiciones entre el orden-desorden de la materia. 

Por otro lado, Giland [8] desarrolla un modelo matemático para un par de ingenieros 
de ecosistemas que se encuentran comúnmente en las tierras secas los cuales forman 
patrones de vegetación y cianobacterias que forman costras del suelo, el cual es un enfoque 
distinto al de HilleRisLambers [12], quien revela que la formación de patrones puede ocurrir 
en áreas semiáridas dada solo la retroalimentación positiva entre la densidad de plantas 
y la infiltración de agua local, junto con la redistribución espacial del agua de escorrentia. 
Estos modelos son algunos de la gran cantidad de hipótesis que se han sugerido sobre el 
origen de formación de patrones de vegetación; no obstante, aún no se ha alcanzado un 


consenso generalizado. Inclusive, se han propuesto modelos estocásticos antagónicos a 
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los modelos de reacción-difusión [4, 26]. 

El objetivo de este trabajo es estudiar analítica y numéricamente un modelo 
matemático propuesto a partir de los modelos de Hardenberg y Lefever, el cual tendrá como 
objetivo describir la distribución vegetal de biomasa, dando así una forma de entender la 
formación de patrones de vegetación y desertificación. Por ello, será necesario considerar 
algunos factores ambientales que tengan causas capaces de explicar la formación de 
dichos patrones. Existe una gran complejidad en las comunidades vegetales, sumado a 
las escalas de tiempo muy largas que rigen la evolución de la ecologia vegetal (respecto 
a los estândares humanos), idear modelos capaces de determinar qué efectos tiene una 
multitud de factores ambientales en la distribución de la vegetación, es una tarea altamente 
problemática y es aún más problemático cuantificar estos efectos. 

El modelo de dinâmica de la distribución vegetal presentado en este documento 
pretende Ilenar los vacíos teóricos. Permitirá predecir que las distribuciones de vegetación 
pueden sufrir una inestabilidad que rompe la simetria como resultado de las interacciones 
entre los diferentes tipos de biomasa, haciendo notar que corresponde al âmbito de los 
fenómenos no lineales y la hidroecologia. Los parâmetros adecuados para describir esta 


inestabilidad corresponden más a las propiedades de la vegetación que al medio ambiente. 


21 EL MODELO DE HARDENBERG - LEFEVER 


Se propone un modelo que emplea ecuaciones diferenciales parciales y no lineales. 
La primera ecuación de este sistema modela propiedades de n(f,t) como la forma de 
su reproducción, la difusión espacial de la misma y las fuentes de recursos que poseen 
dependiendo el entorno en el que se encuentran. Además, la segunda ecuación del modelo 
considera las propiedades que tiene w(7,t) tales como, la evapotranspiración de las 
plantas, una fuente de saturación de agua y la difusión espacial de «w(T,t) influyendo 
sobre n(7,t), Esta variable representa la densidad de la distribución de biomasa y w(7,t) 
representa la humedad, dándole así al modelo un enfoque no lineal más completo. La 


primera de las ecuaciones del modelo se escribe como: 


D=M-Kn+gn?-nº +[Vên+ÃPV*n (1) 


donde el término M está asociado a una fuente/perdida de densidad de biomasa 
estacionaria, el término Kn explica la mortalidad por los herbívoros, nº impide la separación 
de fases, [V*n es la difusión espacial de biomasa y el término IV'm representa un 
término de difusión de largo alcance de la densidad de biomasa [17]. 

Como se busca un modelo más realista, consideramos que la cantidad de agua 
que tienen las plantas no es constante, detalle que nos lleva a acoplar la ecuación para la 
humedad al modelo de Hardenberg [11] a la ecuación (1) de la densidad de biomasa del 


modelo de Lefever [17], de manera que: 
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an 
dt 


= =p=-(1- pow — ewln + ôVÊ(w — Bm) = vo (w— an) (3) 


=-M-kn+gn?-nº + ewin+TVn+=DV'm (2) 


donde p es la precipitación anual; -(1-pn)w es el término de pérdida por evaporación, 
6 V(w-Bn) es el término que modela la difusión de agua mediante absorción de raíces de 
la planta, ew*r) representa el término de acoplamiento que muestra E ganancia/pérdida 
por transpiración del agua en las plantas. Finalmente, el término vp (a) — (tn) muestra 
la pérdida de agua por filtración debido a la pendiente del suelo v, la cual puede ser una 
función espacial [11]. 


31 RESULTADOS 


A continuación, se muestran los resultados numéricos obtenidos mediante el 
método de diferencia finita y aproximando las derivadas espaciales de segundo orden con 
derivada numérica centrada [6, 25, 29], empleado para resolver las ecuaciones (2) y (3) del 
modelo de Hardenberg — Lefever. Se realizaron experimentos numéricos con arreglos de 
mallas de 256x256 a 100,000 iteraciones con un paso de tiempo dt = 0.001, variando los 
parâmetros de control K con un paso dK=0.5 y e con un paso de=0.5. La implementación 
del método numérico fue realizada en lenguaje de programación Python en su versión 3.9 
y sus librerías (Tkinter y numpy) [16]. 

Se encontraron 6 fases fundamentales en las cuales, se puede observar que la 
variable m del modelo significa simultâneamente densidad de biomasa para valores positivos 
y desertificación para valores negativos. El parâmetro de control K del modelo propuesto 
indica la cooperatividad del sistema, mientras que el parâmetro de acoplamiento e modula 
los términos tipo Lotka-Volterra. Los valores de los parâmetros de control restantes son 
fijos y tales que: M=0.1, p=0.4, g=0, [=0.88, y=0 y p=2. Es importante mencionar que, al 
cambiar el valor de M, p y p los resultados que se obtienen no cambian significativamente, 
mientras que, al cambiar los valores de g y y las soluciones numéricas del modelo de 
Hardenberg — Lefever cambian drásticamente, llevando a este a una fase de desertificación. 
Las combinaciones explícitas de los valores de los parâmetros de control, empleados para 
obtener las soluciones numéricas en forma de patrones de vegetación-desertificación, se 
presentan en las Tablas 1 y 2 del Apêndice A. 

Por observación directa de los patrones de vegetación y desertificación, permite 
confirmar las fases de distribución de biomasa — desertificación del modelo Hardenberg — 
Lefever, ya que, además de contar con las fases de puntos (ver Figura 3), huecos (ver Figura 
5), laberintos (ver Figura 4) y fase de vegetación uniforme (ver Figura 6) reportados en las 
referencias [15-28, 7, 11,13, 17, 19, 27, 3-20], además se obtiene una fase entremezclada 
de desertificación con laberintos de mayor complejidad morfológica (ver Figura 2). 

Como una aplicación y resultado final, se realiza una comparación de los resultados 
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obtenidos por nuestro modelo con fotografias capturadas de la formación de zonas de 
biomasa — desertificación en distintas regiones, en las cuales, podemos observar algunas 
de las fases presentadas en las Figuras 1 a 10, tales como, puntos, huecos, laberintos y 
algunas de estas fases son entremezcladas. 

En todos los resultados obtenidos de la simulación computacional, consideramos 
que la pérdida de agua por filtración debido a la pendiente del suelo era nula, es decir, que 
v=0, de manera que, el término vÃ(w — «n) en la ecuación (5) no había sido tomado en 
cuenta. Por esta razón, realizamos un experimento numérico más, en donde consideramos 
los valores de los parâmetros para los arreglos de mallas 256x256 de las fases de puntos, 
laberintos, huecos y la fase entremezclada, con la diferencia de que en este experimento se 
encendió el valor de v arrojando como resultado lo que se reporta a continuación. 





Fig. 1 Fase de desertificación. a)Patrón de distribución de biomasa, b) Patrón de distribución de 
humedad. c)Panorama real ubicado en 3º13'38"N 75º10'059W. 


Fuente: Google Earth. 
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Fig. 2 Fase intermezclada de desertificación con laberintos. a)Patrón de distribución de biomasa, b) 
Patrón de distribución de humedad. c)Panorama real ubicado en 13º20'55"N 2º36'39"E. 


Fuente: Google Earth. 
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Fig. 3 Fase de puntos. a)Patrón de distribución de biomasa, b) Patrón de distribución de humedad. c) 
Panorama real ubicado en 12º57'06"N 3º21'08"E. 


Fuente: Google Earth. 








Fig. 4 Fase de laberintos. a)Patrón de distribución de biomasa, b) Patrón de distribución de humedad. 
c)Panorama real ubicado en 12º38'59"N 3º10'09"E. 


Fuente: Google Earth. 





Fig. 5 Fase de huecos. a)Patrón de distribución de biomasa, b) Patrón de distribución de humedad. c) 
Panorama real ubicado en 19º01'59"S 13º20'28"E. 


Fuente: Google Earth. 
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Fig. 6 Fase de vegetación uniforme.a)Patrón de distribución de biomasa, b) Patrón de distribución de 
humedad. c) Panorama real ubicado en 3º23'31"S 62º10'23'W. 


Fuente: Google Earth. 
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Fig.7 Fase intermezclada de desertificación con laberintos. a)Patrón de distribución de biomasa, b) 
Patrón de distribución de humedad. c) Panorama real ubicado en 13º10'30"N 2º34'52"E. 


Fuente: Google Earth. 
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Fig.8 Fase de puntos. a)Patrón de distribución de biomasa, b) Patrón de distribución de humedad. c) 
Panorama real ubicado en 13º01'36"N 2º06'05"E. 


Fuente: Google Earth. 





Fig.9 Fase de laberintos. a)Patrón de distribución de biomasa, b) Patrón de distribución de humedad. 
c) Panorama real ubicado en 13º07'02"N 2º3816"E. 


Fuente: Google Earth. 





Fig.10 Fase de huecos. a)Patrón de distribución de biomasa, b) Patrón de distribución de humedad. c) 
Panorama real ubicado en 3º12'41"N 75º10'04'W. 


Fuente: Google Eath. 
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41 DISCUSIÓN 


El modelo de Hardenberg-Lefever propuesto, permite entender la formación de 
patrones de vegetación de una manera distinta a la mencionada en la diversidad de modelos 
matemáticos presentados hasta la actualidad. Este modelo es más general y más allá que 
los sistemas de reacción-difusión, cuyas características observables son la densidad de la 
distribución de biomasa y aridez (n) y la humedad (mw). 

La mayoria de las fases mencionadas anteriormente se han observado en diversos 
tipos de ecosistemas reportados en la literatura [1, 7, 17, 20]. Se han encontrado patrones 
de puntos, rayas y bandas en una amplia gama de áreas geográficas en todo el mundo [1, 
7, 17, 20]. Muy a menudo, los patrones involucran a más de una especie. La irregularidad 
autoorganizada en estos ecosistemas difiere en escala y forma. Un ejemplo de patrones 
que han sido reportados en la actualidad son huecos, laberintos, rayas (“arbusto de tigre”) 
y manchas (“arbusto de leopardo”) [24]. 

El mecanismo general que subyace a esta irregularidad autoorganizada es una 
retroalimentación positiva entre el crecimiento de las plantas y la disponibilidad de agua. 
Una mayor densidad de vegetación permite una mayor infiltración de agua en el suelo 
(debido a la penetración de las raíces) y una menor evaporación del suelo (debido a 
la sombra). Como resultado, la vegetación persiste una vez presente. Sin embargo, el 
suelo desnudo es demasiado hostil para la recolonización después de que la vegetación 
desaparezca. En términos generales, esto se debe a inestabilidades inducidas por difusión 
en el mejor de los casos [14, 22], aunque existen otros tipos de inestabilidades que ocurren 
mediante otros mecanismos de inestabilidad [2, 14, 21, 23]. 

La formación de patrones es parte de la dinâmica inestable de las ecuaciones 
no lineales dadas por el modelo Hardenberg — Lefever, como producto de una dinâmica 
inestable. La solución numérica del modelo nos muestra una diversidad de diferentes 
fases morfológicas en los patrones. Esta inestabilidad es causada por la competencia de 
distribución de biomasa sobre los recursos hídricos. El tipo de patrón depende del rango 
de acumulación de humedad por planta; fase uniforme de vegetación para valores de £ 
relativamente altos, laberintos (o franjas/tripas) para cuando el valor de £ toma un valor 
intermedio y agujeros en una cobertura uniforme para valores de £ relativamente bajos. 
Estos patrones reflejan una auto-organización óptima del sistema con respecto a los 
recursos hídricos, esto se sabe debido a la dinâmica no lineal y la competencia interacción 
de los términos del modelo. 

De acuerdo con Mander y Dekker [19], estos patrones pueden observarse en 
imágenes satelitales y pueden producirse mediante sistemas de activación-inhibición en 
modelos computacionales. O bien, mediante software que permite una caracterización de la 
distribución de biomasa o desertificación de estas imágenes satelitales [5], lo cual es motivo 
de estudios futuros. El desarrollo de patrones de vegetación espacial en simulaciones sigue 
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una secuencia bien establecida que está relacionada con la cantidad de Iluvia suministrada 
a la superficie terrestre. Con niveles de Iluvia relativamente altos, toda la superficie de la 
tierra está cubierta de vegetación y, a medida que la Iluvia disminuye progresivamente, 
los patrones de vegetación cambian de huecos (cubierta de vegetación casi continua con 
pequenas aberturas) a laberintos (redes reticuladas de vegetación) a manchas (pequerfos 
parches de vegetación), y finalmente al suelo desnudo, tal y como se puede observar en 
las Figuras 1 - 6. 

El aumento de la infiltración de agua debajo de los patrones de vegetación y la 
sombra representan importantes procesos de facilitación local, mientras que el flujo de 
agua terrestre y la competencia a través del sistema de raíces son ejemplos de importantes 
retroalimentaciones a largo plazo. Kletter y von Hardenberg [15], mencionan que, un 
resultado importante de estos estudios es su robustez y consistencia: todos los modelos 
son capaces de reproducir la misma fenomenologia de patrones periódicos de biomasa tan 
pronto como se eligen rangos realistas para los parâmetros. Se diferencian principalmente 
en el detalle con el que representan diferentes mecanismos de retroalimentación, lo que 
permite a algunos de ellos capturar con mayor detalle la fenomenologia de los patrones de 
vegetación, su resiliencia a los cambios ambientales o la creación de nichos ecológicos. 

Presentamos aquí un nuevo modelo matemático que da cuenta de una amplia gama 
de patrones de biomasa — desertificación. Usando este modelo demostramos la ventaja 
potencial de un enfoque de la dinâmica no lineal espacio-temporal en el estudio de los 
ecosistemas (Ecohidrologia). La coexistencia de estados estables e inestables arroja 
nueva luz sobre los fenómenos de desertificación y vegetación. La secuencia de estados 
estables a lo largo de un gradiente de Iluvia motiva una clasificación de la aridez que refleja 
mejor las propiedades del sistema, como la vulnerabilidad a la desertificación. 


5 | CONCLUSIONES 


En general, estos patrones de vegetación se interpretan como un proceso de 
autoorganización descrito por mecanismos no lineales dados por las ecuaciones del modelo 
Hardenberg — Lefever. La dinámica de separación de fases encontrada y sus transiciones 
entre sus fases, refleja el cambio de acuerdo con los parâmetros sensibles e y K de la 
estructura de la distribución de biomasa frente al estrés hídrico. 

El modelo de Hardenberg — Lefever es un modelo bastante completo y competente, 
pues además de mostrar en los resultados una cantidad considerable de reproducción de 
patrones de vegetación y desertificación de diferentes ecosistemas, también puede ser 
comparado con la naturaleza de manera inmediata mediante imágenes satelitales. Mediante 
ésta comparación de resultados, se ha demostrado que el modelo de Hardenberg — Lefever 
es un modelo descriptivo. Estas observaciones muestran que pueden surgir patrones auto- 
organizados a diferentes escalas [7, 8, 11. 18-24]. 
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A diferencia de otros modelos existentes, tales como el modelo de Rietkerk [24], el 
modelo de Lefever — Lejeune [18] o el modelo de Siero [27], los cuales utilizan sistemas de 
ecuaciones diferenciales parciales más complicados, el modelo de Hardenberg-Lefever es 
más fácil de tratar analítica y numéricamente. 

En el contexto de los patrones de biomasa — desertificación espacial en ecosistemas 
de tierras secas, el enfoque que hemos desarrollado, se suma a la creciente caja de 
herramientas de métodos para analizar estos patrones. En este contexto, el trabajo futuro 
podríia explorar el grado en que nuestra caracterización de la morfologia del patrón de 
biomasa podria usarse como una herramienta de detección temprana del estado de 
desertificación. 

Finalmente se demuestra que el patrón de vegetación es conceptualmente una 
transición desorden — orden, cuya dinâmica es similar a los procesos no lineales de 
separación de fases, más que a la inestabilidad inducida por difusión. 
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APÊNDICE A: TABLAS CON LOS VALORES DE LOS PARÁMETROS DE 
CONTROL PARA LAS SOLUCIONES NUMERICAS DEL MODELO DE 
LEFEVER-HARDERG 


Parameters Fig. 1(a,b) Fig. 2(a,b) Fig. 3(a,b) Fig. 4(a,b) Fig. 5(a,b) Fig. 6(a,b) 


K 0.52 1.34 0.45 0.45 0.44 1.14 
£ -0.67 0.115 -0.67 -0.88 1.28 -0.42 
p 2 2 2 2 2 2 
M 0.1 0.1 0.1 0.1 0.1 0.1 
p 0.4 0.4 0.4 0.4 0.4 0.4 
r 0.88 0.88 0.88 0.88 0.88 0.88 
g 0 30 0 0 0 0 

Grid 256º 256º 256º 256º 256º 256º 


Tabla 1: valores de los parâmetros de control de las figuras 1 a 6. 
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Parameters Fig. 7(a,b) Fig. 8(a,b) Fig. 9(a,b) Fig. 10(a,b) 

K 1.34 0.45 0.45 0.44 
E 0.115 -0.67 -0.88 1.23 
p 2 2 2 2 
M 0.1 0.1 0.1 0.1 
p 0.4 0.4 0.4 0.4 
r 0.88 0.88 0.88 0.88 
g 30 0 0 0 

Grid 256º 256º 256º 256º 


Tabla 2: valores de los parâmetros de control de las figuras 7 a 10. 
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RESUMO: As formigas cortadeiras são 
responsáveis por prejuízos agrícolas. Seu 
controle é basicamente químico, o que causa 
prejuízos ambientais, aos animais selvagens 
e a saúde do homem. Medidas preventivas de 
controle devem ser tomadas, e o uso de plantas 
inseticidas é uma alternativa viável. Desta 
maneira, o objetivo deste trabalho foi investigar o 
potencial de ação inseticida das folhas de laranja 
(Citrus sinensis L.) e Bt (Bacillus thuringiensis) 
sobre as formigas cortadeiras (A. sexdens 
sexdens) em condições experimentais. O projeto 
foi executado no laboratório de zoologia do 
Instituto de Saúde e Biotecnologia — ISB/UFAM, 
Campus/Coari. A coleta das formigas cortadeiras 
foi realizada no Centro de Apoio à Pesquisa do 
Médio Solimões, localizado na cidade de Coari. 
As folhas de laranja foram coletadas no campus 
do ISB, encaminhadas ao laboratório de química, 
onde foi realizada assepsia com água destilada. 
As amostras foram separadas, pesadas e 
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levadas à estufa de circulação de ar forçada a 
uma temperatura de 55º C, por um período de 
cinco dias. As folhas foram trituradas e levadas 
ao sistema Soxhlet, tendo como solvente água 
destilada. No bioensaio, foram colocadas 10 
formigas adultas para cada placa de Petri, com 
auxílio de pinças entomológicas. Posteriormente, 
foi aplicado cerca de 3 mL dos extratos botânicos 
em cinco concentrações (1:0; 1:42; 11; 114%; 
1:2). A testemunha foi pulverizada com água 
destilada. Foram utilizados cinco tratamentos com 
três repetições, mais a testemunha totalizando 
18 parcelas para cada produto. O experimento 
foi observado diariamente durante cinco dias. 
Os dados foram submetidos a uma análise de 
variância com 95% de significância. Registramos 
95% de mortalidade das formigas nas 120 horas 
de observação em todos os tratamentos, mas 
não houve diferença significativa. Novos testes 
estão sendo realizados para verificar a eficiência 
dos produtos vegetais. 

PALAVRAS-CHAVE: Controle 
plantas inseticidas, insetos pragas. 


alternativo, 


USE OF BOTANICAL EXTRACTS OF 
ORANGE LEAVES (Citrus sinensis 
L), ON HARVESTING ANTS (Atta 
sexdens sexdens) IN EXPERIMENTAL 
CONDITIONS 


ABSTRACT: Leaf-cutting ants are responsible 
for agricultural losses. Its control is basically 
chemical, which causes damage to the 
environment, wild animals and human health. 
Preventive control measures must be taken, and 
the use of insecticidal plants is a viable alternative. 
Thus, the objective of this work was to investigate 
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the insecticidal action potential of orange (Citrus sinensis L.) and Bt (Bacillus thuringiensis) 
leaves on leaf-cutting ants (A. sexdens sexdens) under experimental conditions. The project 
was carried out in the Zoology laboratory of the Institute of Health and Biotechnology — ISB/ 
UFAM, Campus/Coari. The collection of leaf-cuiting ants was carried out at the Middle 
Solimões Research Support Center, located in the city of Coari. The orange leaves were 
collected on the ISB campus, sent to the chemistry laboratory, where asepsis was performed 
with distilled water. The samples were separated, weighed and taken to a forced air circulation 
oven at a temperature of 55ºC for a period of five days. The leaves were crushed and taken to 
the Soxhlet system, using distilled water as solvent. In the bioassay, 10 adult ants were placed 
in each Petri dish, with the aid of entomological tweezers. Subsequently, approximately 3 mL 
of botanical extracts were applied in five concentrations (1:0; 1:72; 1:1; 1:1%; 1:2). The control 
was sprayed with distilled water. Five treatments were used with three replications, plus the 
control, totaling 18 plots for each product. The experiment was observed daily for five days. 
The data were submitted to an analysis of variance with 95% of significance. We recorded 
95% ant mortality in the 120 hours of observation in all treatments, but there was no significant 
difference. New tests are being carried out to verify the efficiency of plant products. 
KEYWORDS: Alternative control, insecticidal plants, insect pest. 


INTRODUÇÃO 


As formigas cortadeiras são responsáveis por grandes prejuízos econômicos, sendo 
que as formigas pertencentes aos gêneros Aita (saúvas) e Acromyrmex (quenquéns), 
são as que mais se destacam como um dos principais insetos-pragas, que atacam os 
sistemas agrícolas, nas regiões de clima tropical e subtropical das Américas. Dentre estas, 
as sauúvas são as maiores causadoras de danos, devido ao grande tamanho populacional 
das colônias, desta forma para suprir a necessidade de alimentação destes insetos, as 
formigas precisam de grandes quantidades de folhas, que servem de substrato para o 
fungo simbionte que compõe a principal fonte de alimentação do sauveiro (HOWARD et 
al., 1988) 

Existe uma forte relação entre as formigas cortadeiras e o fungo simbionte, pois a 
sobrevivência de ambos está relacionada a essa associação, e acredita-se que estes não 
podem sobreviver separadamente. As formigas são responsáveis pelo fornecimento de 
substrato para o fungo, além de manterem o ambiente limpo e livre de contaminantes, o 
que favorece o desenvolvimento do fungo (HOLLDOBLER; WILSON, 1990). As formigas 
cortadeiras são de grande importância econômica pelo fato de utilizarem grande quantidade 
de material vegetal para o cultivo do fungo simbionte. 

Segundo Della Lucia e Oliveira (1993), o processo de obtenção de material vegetal 
para suprir as necessidades alimentar da colônia, possui alguns critérios como, seleção do 
material, corte e transporte das folhas para interior do ninho. Para a eficiência das atividades 
de forrageamento, as formigas possuem estratégias comportamentais para encontrar 


e utilizar fontes energéticas e nutrientes, isso resulta na interação de comportamentos 
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individuais das operárias e na construção de trilhas de exploração (SCHLINDWEIN, 2004). 

Durante as atividades de forrageamento as saúvas, exploram enormes áreas, 
fazendo o corte em diversas espécies de plantas, ao longo dessas áreas são construídas 
trilhas, de aproximadamente 400 metros de comprimento por 20 cm de largura, essa trilha 
é usada para o transporte do material vegetal. As formigas possuem hábitos noturnos para 
as atividades de forrageamento, porem essas atividades podem ser observadas durante 
o dia no período da manhã com o clima mais ameno. Apesar de se alimentarem de varias 
espécies vegetais, as formigas apresentam características bastante seletivas, possuem 
preferencias com relação ao corte das plantas, por partes mais jovens e macias, como 
folhas novas, brotos e flores (LIMA et al., 2001). 

As formigas são conhecidas pela sua voracidade, em relação ao ataque que fazem 
nos sistemas agroflorestal, os danos causados por essas formigas variam de danos parciais 
a total desfolhamento das plantas, as ações das formigas são intensa e constante podendo 
acontecer em qualquer fase de desenvolvimento das culturas. Por estes motivos as 
formigas cortadeiras são denominadas uma das principais pragas das culturas (ZANETTI 
et al., 2002). 

Estes insetos representam uma forte ameaça às plantações, sendo responsável 
por sérios danos econômicos a agricultura nacional, pelo fato da competição pelas plantas 
cultivadas pelo homem. As formigas cortadeiras estão distribuídas por quase todo o 
território brasileiro, e os ataques são realizados intensamente durante o ano inteiro. Outra 
questão preocupante está relacionado às dificuldades e os altos custo no controle destas 
formigas. Para muitos agricultores as formigas cortadeiras estão na linha de frente como os 
principais inimigos das culturas, e para o efetivo controle, é necessário muito investimento 
com formicidas, equipamentos e mãos de obra (LOECK; GRUTZMACHER, 2001). 

Em virtude dos danos econômicos causados pelas formigas cortadeiras, diversas 
formas de controle já foram desenvolvidas (JACCOUD, 2000). Apesar dos vários métodos 
de controle existente, os métodos mais utilizados no controle das formigas cortadeiras são 
uso de inseticidas sintéticos, porem estes podem gerar impactos negativos ao ambiente e ao 
homem (SOUZA-SILVA et al., 2005). Segundo Boaretto e Forti (1997) o controle químico de 
formigas cortadeiras é realizado por meio de iscas, termonebulização e fumigantes, contudo 
estes métodos exigem altos custos. Desta forma grandes quantidades de inseticidas são 
aplicadas na tentativa de controle dessas formigas e, nem sempre apresentam resultados 
satisfatórios (DELLA LUCIA; FOWLER, 1993). 

Para Júnior (2003) é de fundamental importância à busca por novas alternativas, que 
posam fazer frente ao controle de pragas, e que o mesmo possa oferecer maior segurança, 
seletividade, biodegradabilidade, viabilidade econômica e aplicabilidade e baixo impacto 
ambiental em programas integrados de controle de insetos. Os diversos efeitos negativos 
advindo do uso indiscriminado de agrotóxicos sintéticos podem ser minimizados com o 


uso de produtos alternativos como os extratos botânicos, devido a algumas características 
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benéficas relativas à toxicidade, pois possuem degradação rápida o que geralmente pode 
reduzir o impacto aos inimigos naturais, ao homem e ao ambiente (BRITO et al., 2008). 

Muitos estudos estão sendo desenvolvidos sobre plantas com potencial inseticida, 
com o objetivo de extrair e identificar compostos que possam atuar no controle das 
formigas ou a seu fungo, esses estudos têm sido realizados com o intuito de desenvolver 
metodologias viáveis para aplicação em campo. O uso de substâncias tóxicas, extraídas 
de vegetais, poderá apresentar-se, no futuro, como uma possibilidade de controle (ACÁCIO 
BIGl et al., 2004). Logo esta pesquisa teve como objetivo investigar o potencial de ação 
inseticida das folhas de laranja (Citrus sinensis L.) sobre as formigas cortadeiras (A. 
sexdens sexdens) em condições experimentais e comparar a toxicidade da mortalidade 
das concentrações dos extratos das folhas de laranja. 


MATERIAIS E MÉTODOS 


O experimento (bioensaios) foi realizado nos laboratórios de Química Analítica, 
Química Geral e Química Orgânica, e zoologia localizados na Universidade Federal do 
Amazonas na cidade Coari. 


Coleta dos insetos 


A coleta das formigas cortadeiras foi realizada no Campus UFAM / ISB, localizados 
na cidade de Coari, nas coordenadas S 04º05'37.5 W 063º09'04.1. 
Coleta da planta 


Amostras de folhas de laranja (Citrus sinensis) foram coletadas no campus do ISB/ 
UFAM, encaminhadas ao laboratório de química, no qual foi realizada uma assepsia com 
àgua destilada. As folhas prensadas para a confecção da exsicata com fins de identificação 
e catalogação e encaminhadas ao laboratório de botânica do ISB. 


Preparo das concentrações 


No laboratório de Química Orgânica as folhas dos vegetais citados acima foram 
separadas, pesadas e levadas a estufa de circulação de ar forçado a uma temperatura de 
55ºC, por um período de cinco dias. As folhas retiradas da estufa, já secas (500 9) foram 
trituradas. O material triturado foi encaminhado ao sistema Soxhlet, tendo como solvente 
àgua destilada. 


Bioensaio: exposição via de contato 


Foi realizado o ensaio com objetivo de testar a eficiência do método de aplicação do 
produto: exposição via de contato. 

Foram colocados 10 formigas adultas, com idade de 10 a 15 dias, para cada placa 
de Petri, com auxílio de pinças entomológicas. Após foi aplicado cerca de 3 mL do extrato 
botânico nas concentrações de (1:0; 1/1/2; 1:1; 1:1 Y%; 1:2), totalizando 5 concentrações 
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(cinco tratamentos) com cinco repetições, mais a testemunha, totalizando 30 parcelas. As 
observações foram diárias, totalizando 120 horas ou cinco dias consecutivos de visualização 


do experimento. A testemunha foi pulverizada com água destilada. 


Estatística e análise de dados 


Os dados foram submetidos a uma análise de variância — ANOVA (Zar, 1984), 
onde foi usado um delineamento inteiramente casualizado com cinco tratamentos e cinco 


repetições, sendo avaliada a mortalidade dos insetos. 


RESULTADOS E DISCUSSÃO 


As observações referentes à taxa de mortalidade foram realizadas a cada 06 horas 
durante 5 dias (120 horas). Os valores percentuais referentes à mortalidade dos insetos, 
nas concentrações de (1:0;1/1/2; 1:1; 1:1 4; 1:2) do extrato de laranja e testemunha, no 


intervalo de tempo de 120 horas, encontra-se no gráfico 01. 
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Gráfico 01 - Mortalidade das formigas cortadeiras, nas 120horas de observação. 


De acordo com GALLO et al. (2002), o objetivo principal do uso de extratos vegetais 
é reduzir o crescimento da população de pragas. Segundo os autores, a mortalidade do 
inseto é apenas um dos efeitos e que, geralmente, necessita de concentrações muito 
elevadas. 

Segundo VENDRAMIM (1997) as avaliações podem ser feitas em campo, casa 
de vegetação ou em condições de laboratório. Em campo ou em casa de vegetação, 
utilizando-se parcelas tratadas ou não, o efeito pode ser determinado através da avaliação 
da população e oviposição do inseto ou do dano sofrido pela planta. Em condições de 
laboratório, quando são oferecidas folhas provenientes de plantas tratadas e não tratadas, 
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são avaliados a oviposição, consumo de alimento (em testes com e sem chance de 
escolha), duração do ciclo biológico, peso e tamanho, mortalidade das fases imaturas e da 
fase adulta, fecundidade, fertilidade e alterações morfogenéticas. 

O uso de extrato vegetais ressurge como uma opção para o manejo integrado 
de pragas onde associado a outras práticas, pode contribuir para a redução de doses 
e aplicações de inseticidas químicos sintéticos, além de amenizar a mortalidade dos 
organismos benéficos e ao meio ambiente. 

Os testes com os extratos de laranja apresentaram uma mortalidade de 100% dos 
indivíduos nas 120 horas de observação. Na concentração 1:0 (extrato puto) após 24 horas 
dos testes 100% das formigas estavam mortas nas cinco repetições. Na concentração 
de 1/2 a taxa de 100% de mortalidade dos insetos foi alcançada após 48 horas do início 
dos testes. As demais concentrações (1:1; 1:1 Y; 1:2) após 72 horas possuíam 100% dos 
insetos mortos em todas as placas. Não ocorreu mortalidade da testemunha em todo o 
experimento. 

Lima et al (2009) testou atividade inseticida do óleo essencial de Piper hispidinervum 
sobre lagarta-do cartucho do milho (Spodoptera frugiperda) por ingestão e contato tópico. 
No Brasil a lagarta-do-cartucho Spodoptera frugiperda (J. E. Smith, 1797) (Lepidoptera: 
Noctuidae), é considerada a principal praga da cultura de milho, atacando plantas jovens e 
reduzindo em até 34% a produção; dependendo da idade da planta ataca também a cultura 
da cana-de-açúcar, o arroz e o algodoeiro (VALICENTE & CRUZ, 1991). Os bioensaios 
demonstraram sua atividade inseticida para S. frugiperda, causando mortalidade e redução 
alimentar pelo teste de ingestão e toxicidade aguda pelo teste de aplicação tópica, sendo 
também observados sintomas de neurotoxicidade, como agitação e hiperatividade. 


CONSIDERAÇÕES FINAIS 


Na concentração 1:0 (extrato puro) ocorreu a mortalidade de 100% dos indivíduos 
nas primeiras 24h. Todas as concentrações utilizadas mostraram-se eficientes no controle 
das formigas cortadeiras em condições experimentais. No entanto, novos testes precisam 
ser realizados em campo, para recomendar tal produto para agricultores. 

Avanços nessa linha de conhecimento é de fundamental importância para auxiliar no 


controle das pragas para produtores que são adeptos à agricultura orgânica (agroecologia). 
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